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mri;TC*V'-r.  *J  jxk  A TURBULENT  TO  A LAMINAR  BOUNDARY  LAYER 


ABSTRACT 

The  -Jundary  layer  on  a $&'J  included  angle  cone -cylinder  model  has  been 
studied  ever  * range  of  supersonic  npeods . Boundary  layer  recovery  factors 
have  been  measured  and  shadowgraphs  have  been  taken  of  the  model  in  the  wind 
tunnel  and  free  flight.  These  measurements  indicate  that  when  the  boundary 
layer  ahead  of  the  shoulder  is  turbulent,  the  boundary  layer  is  laminar  after 
the  shou]  ier  for  a short  distance  . A theoretical  model  of  the  phenomenon  has 
bean  constructed  and  calculations  have  been  made  using  available  boundary 
laysr  t-hsory  amd  4 vpayj  ih/w>ta.*t  data  to  describe  the  course  of  t.h.*?  bo'-wtay 
layer  frem  .lust  ahead  of  the  shoulder  to  the  transition  back  to  turbulent  flow 
on  the  cylinder.  xteascnaoiy  gooci  agreement  is  fourn  U G vWyib*i  tiiS  uXp6r XXQB Iiw5 
and  the  calculations* 


A distinction  i3  made  between  the  complete  shear  layer  at  a given  position 
which  results  from  the  action  of  frictional  forces  all  along  the  model  surface 
frail  the  model  tip,  and  the  boundary  layer  at  the  same  position  which  may  be  all 
or  only  the  inner  portion  of  the  complete  shear’  layer, 

A turbulence  parameter  analogous  to  Taylor* s turbulence  parameter  fee  the 
effect  of  a turbulent  free  stream  on  laminar  rrbulent  transition,  but  appli- 
cable to  a high  speed  boundary  layer  is  obtained  • 
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Recovery  factor 


Laminar  recovery  factor 
Turbulent  recovery  factor 
Wall  temperature 
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to  conditions  just  ahead  of  the  corner 
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Boundary  layer  thickness 
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* Local  friction  coefficient 


• 2G,  - Average  frirt : coefficient 


yT 


frv 

■ Friction  velocity  «,/  — - 

v A.r 


S(M) 


- CF^ 


fW 

v 


i i > 

V.  V.  V 


* coot  mean  square  turbulent  fluctuation  velocities  in  the 
, j . : '?  'erections 


g tv  n 

« -*  - 


■ Initial  coordinates  in  the  Lagraagian  description 
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It  has  been  generally  believed  that  ones  a 'boundary  Isycr  bsccsoesi  tur- 
bulent it  remains  turbulent  * In  tbs  usual  transition  regies,  say  of  a 


boundary  layer  on  a xl»v  plat**,  an  instantaneous  picture  wight  show  laminar 
fact  free  the  leading  edgo,  a generally  laminar  region  interspersed  with 


"bursts"  of  turbulent  fluid,  »r>d  then  finally,  further  downstream  a fully 
developed  turbulent  boundary  layer.  In  this  case,  tha  turbulent  "bursts" 
remain  turbulent  to  they  move  downs  Lrwam  and  grow,  enveloping  the  contiguous 
laminar  regions.  The  fluid  in  the  laminar  patches  between  the  "bursts"  has 
been  laminar  from  the  leading  edge.  In  this  report,  evidence  will  be  pre- 
sented which  indicates  that  under  suitable  conditions  a fully  developed 


turbulent  boundary  layer  can  be  converted  into  a l aminar  boundary  layer , 

The  transf arv.tf  does  not  require  the  disappearance  of  the  turbul  ent  energy 

in  tV  i boundary  layer.  Rather,  same  of  the  turbulent  fluid  in  the 

trubulent  boundary  layer  becomes  laminar  in  the  saw  sense  that  some  of  the 
turbulent  fluid  in  a wind  tunnel  airat^v-i  beewnan  laminar  when  there  is  a 
laminar  boundary  layer  on  the  surface  of  a wind  tunnel  model* 


The  phenomenon  of  inverse  transition  was  noticed  accidentally.  A 58° 

- angle  cone-cylinder  model  was  being  used  to  study  the  effect  of  a 

sharp  proosvu-6  gradient  on  the  boundary  layer  Wapermoure  recovery  factors  * 
The  idea  was  to  produce  reasonably  well-defined  laminar  and  turbulent  bound- 
ary layers  on  the  con?  and  then  measure  the  equilibrium  surface  temperatures 
on  the  cylinder.  We  believed  that  the  velocity  via  temperature  distributions 
in  the  boundary  layers  would  not  be  the  familiar  equilibrium  ones  and  that 
the  surface  temperature  would  be  thereby  changed.  It  was  thought  that  the 
deviation  would  be-  the  greatest  Just  after  the  shoulder  of  the  model  and  that 
the  usual  equilibrium  values  would  be  approached  asymptotically . With  a 
smooth  model  surface,  the  boundary  layer  was  laminar  on.  he  cone  and  became 
turbulent  on  the  cylinder  at  different  positions,  depending  on  the  tunnel 
pressure  level . The  surface  temperatures  on  the  cylinder  were  in  general 
agreement  with  cur  evpnct-t-icss,  rising  in  a normal  way  in  the  transition 
region  on  the  cylinder.  Whan  transition  on  the  cone  was  fixed  by  surface 
roughness,  the  temperatures  on  the  cylinder  followed  a distinctly  different 
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pattern.  For  a short  distance  along  the  cylinder,  the  surface  terape-  'tires 
were  approximately  constant  and  corresponded  closely'  to  the  usual  laminar 
levels  after  which  the  temperatiu’es  rose  rapidly  to  values  close  to  the 
usual  turbulent  levels . 


After  some  unsuccessful  attempts,  5 had  “..graphs  wore  obtained  of  the 
model*  With  surface  roughness  the  fully  developed  turbulent  boundary  layer 
ahead  of  the  shoulder  was  clearly  visible . From  the  shoulder  downstream, 
the  turbulence  that  was  visible  in  the  boundary  layer  ahead  of  the  shcwlder 
seemed  to  disappear,  or  at  least  it  was  not  distinguishable  from  the  back- 
ground of  the  picture.  In  the  downstream  region  vdiere  the  surface  temper- 
atures increased-  a turbulent  boundary  layer  was  again  visible. 


This  report  presents  the  results  or  a study  of  the  effects  of  the  strong 
favorable  pressure  gradient  at  the  shoulder  or  the  cone-cylinder  on  both 
laminar  ard  turbulent  boundary  layers. 


EXPERIMENTAL  ARRANGEMENT 


The  tests  were  conducted  in  the  Ballistic  Research  Laboratories'  Flexible 
Ncssle  Tunnel  shown  in  Figure  1.  The  tunnel  ia  of  the  closed  circuit  variable 
density  type  so  that  data  could  be  obtained  ever  a range  of  Reynolds  Numbers 
at  each  Mach  Number  as  shown  in  the  following  tables 


M 


3e02 

3.55 


\ cm /Max 

2.8  x 10* 

1.8  x 10* 


'Re\ 

^ cm  min 


.70  x 10' 


J;9  x 10' 


me  test  section  at  these  Ma  C**  Nuxfcers  is  38  cm.  wide  by  33  can*  high, 
and  the  top  and  bottom  walls  of  the  test  section  are  sloped  to  account  for 
the  boundary  layer  growth  on  the  tunnel  walls.  The  Mach  Number  and  flow 
direction  distributions  along  the  axis  of  the  tunnel,  including  the  region 
U9*d  for  those  tests,  are  shown  in  Figure  - « 
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for  those  teste  is  shown  in  Figure  3* 
ijurf^Cv  Wo^h  Nuinbs?  dlst-Jibutics  4^  fW«T» 


to  analyse  the  temperature  data,  a steel  model  with  static  pressure  orifices 
along  its  length  (see  Figure  u)  was  tested  at  M * 3*02,  the  Hash  Number  fet 
which  most  of  the  measurements  were  made.  Figure  5 shows  the  conperiaon  of 
the  measured  speed  distribution  along  the  cyD  inder  and  the  surface  Mach 
Wumbers  given  in  reference  1. 


If,  in  calculating  too  recovery  factor,  ar.  incorrect  value  for  the  local 
Mach  Number  is  used,  then  the  computed  recovery  factor  will  be  In  error.  If 
the  Mach.  Number  usan  is  too  high,  then  the  recovery  fact. or  will  be  high  and 
vice  versa  * I;inc«>,  at  M as  3,  * .0 5 wrror  is  Mach  Number  only  restlts  in  an 
error  of  „oui  in  r,  the  sgreeisent  in  Figure  5 was  considered  good  enough  to 
vise  the  calculations  at  the  other  tea;:  Mach  Numbers.  The  calculated  surface 
distribution  for  M = 3.5?  is  shown  in  Figui-e  6. 


The  construction  of  the  Incite  temperature  rec-.  -ery  model  is  shown  in 
Figure  3.  Iron-Cons tantan  thermocouple  leads  -fare  attaeii&d  to  small  lead 
plugs  placed  at  regular  intervals  along  the  lucite  body.  Because  the  Mach 
Number  is  lower  on  the  cone  than  on  the  cylinder,  the  level  of  the  tempera- 
tures on  the  cone  is  higher  than  the  level  on  the  cylinder.  Although  the 
lucite  is  a moderately  good  insulator,  the  model  thickness  was  reduced  near 
the  shoulder  to  reduce  the  heat  transfer  in  the  lucite  from  ahead  of  She 
shoulder  to  the  cylinder.  The  spoked  body  supports,  as  well  as  the  front 
portion  of  the  strut,  vere  also  made  of  lucite.  The  interior  of  the  modal 
was  filled  with  glass  wool  and  the  base  was  covered  with  lucite,  except  for 
a hole  to  take  out.  the  thermocouple  Deads  and  a vent  hole  to  allow  the  in- 
terval pressure  of  the  model  to  equalise  with  model  base  pressure. 


The  shadowgraphs  were  taken  Ly  placing  a spark  source  close  to  the  out- 
side surface  of  one  window  at  the  test  section,  and  the  film  against  the 
other  window.  The  spark  source  was  developed  by  the  Airflow  Branch  of  the 
Exterior  Ballistics  Laboratory  for  taking  free  flight  interferograms,  and 
had  high  intensity  coupled  with  short-  duration.  At  these  high  Mach  Numbers, 
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turbulence  In  the  boundary  layer  sweeps  by  with  almost  olio  %TO&  stress:  vel- 
ocity . Pr^sumabjy,  the  shadowgraph  will  include  the  optical  contribution 
of  the  turbulent  boundary  layers  on  the  windows  as  well  as  the  boundary  layer 
on  the  model,  To  reduce  the  interference  of  the  window  boundary  layers,  the 
light  source  and  the  film  v=>re  planed  close  to  the  outside  of  each  of  the 
te3i  section  windows  suggested  by  Professor  L;  s . G,  Kcwnsenay  of  .Tohns 
Hopkins  University. 
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Probably  the  most  reliable  information  on  idle  boundary  layer  recovery 
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factors  in  the  transition  regj  on  hes  bean  obtained  on  cones 


Figure  7 


shows  some  of  these 


taken  from  t,ve  different  win*’  tunnels  at  various 


Mach  Numbers.  As  can  be  seen  from  the  figure,  the  recovery  factors  rise  from 
about  .81(5  where  the  boundary  layer  is  laminar,  to  turbulent  values  of  about 
.88.  The  curves  shown  cover  a wide  range  of  transition  Reynolds  Numbers, 
presumably  due  to  differences  in  the  disturbance  levels  in  the  various  wind 
tunnels  * As  with  skin  friction  transition  curves,  when  the  transition  starts 
earlier,  tbs  transition  region  on  a logarithmic  3caj.e  tends  to  be  longer. 

SMOOTH  MODEL  RESULTS 


1 I 


For  the  58  cone  cylinder  in  the  smooth  condition,  that  ±3  with  no  tur- 
bulence promoters  on  the  cone,  the  boundary  layer  was  laminar  on  the  cone  and 
for  various  distances  along  the  cylinder.  Figures  8 and  9 are  shadowgraphs 
at  a supply  pressure  of  litO  cmj  the  boundary  layer  is  laminar  on  the  cone  and 
undergoes  transition  on  the  cylinder  at  abexit  20  cm.  from  the  shailder. 
Several  features  of  the  temperature  data  shown  in  Figure  10  are  of  particular 
interest.  The  temperature  level  for  the  first  six  plugs  is  approximately  the 
same,  but  a systematic  decrease  can  be  observed  as  the  shoulder  is  approached 
and  the  tempe^tsT*®  last  chc=d  of  the  shoulder  is  considerably  less  than  the 
others.  On  the  cyliiider,  the  temperatures  apparently  rise  towards  an  asymp- 
totic value  well  back  on  the  cylinder.  The  temperatures  within  one  or  two 
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centimeters  of  the  shoulder  fall  above  the  curve  that-  raj 'osonts  the  moat 
natural  extrapolation  upstream  of  the  temperature  data.  The  temperature  fall 
off  ahead  of  the  shoulder  arid  the  higher  values  just  behind  the  shoulder  are 
believed  to  be  the  xeault  of  heat  transfer  in  the  itself,  since  the 

lucite  necessarily  must  have  a continuous  temperature  distribution.  The  same 
data  are  shown  as  recovery  factors  in  Figure  11.  The  recovery  factors  are 
lower  than  the  normal  laminar  values  after  the  shoulder,  rise  smoothly  towards 
an  asymptotic  value,  until  there  is  an  abrupt  rise  at  the  onset  of  transition 
on  the  cylinder.  At  first  it  was  thought  that  heat  transfer  effects  in  the 
vicinity  of  the  shoulder  might  be  responsible  for  the  low  temperatures  after 
the  shoulder.  Since  the  lass  of  the  cons  is  cooled  by  heat  transfer  to  the 
front  of  the  cylinder,  there  must  be  heat  transfer  from  the  boundary  layer 
ahead  of  tho  shoulder  to  the  cons  surface.  After  the  corner,  the  boundary 
layer  would  iiave  an  energy  deficiency  which  would  result  in  i recovery 
temperatures.  However,  if  it  is  assumed  that  the  heat  flow  can  be  mainly 
accounted  for  by  two  mechanisms } heat  transfer  between  the  model  and  the  air, 
and  heat  flow  within  the  shell,  (heat  flow  to  the  air  within  the  model  is 
neglected),  then,  since  the  model  is  at  equilibrium  temperature  at  all  points, 
heat  flow  ou~.  t-ne  eir  ah:-i  cf  the  shoulder  should  be  balanced  by  heat  flew 
back  to  the  air  sf  U the  shoulder  so  that  the  heat  flow  in  the  Incite  should 
not  cause  an  energy  deficiency  the  boundary  layer  on  the  blinder.  The 
heat  flow  out  of  and  back  into  tha  boundary  layer  would  be  expected  to  have 
some  influence  on  the  temperature  distribution,  and  so  the  recovery  temper- 
atures. Heat  flow  estimates  in  the  lucite  indicate  that  these  heat  transfer 
effects  should  be  negligible  2 .F?  cm.  from  the  shoulder  at  Pq  ■ lUO  cm. 

The  laminar  recovery  factors  on  the  cone  are  in  good  agreement  with  the 
data  shown  in  Figure  7.  In  the  initial  tests,  a great  deal  of  dlffic  - ■ was 
experienced  in  repeating  tha  recovery  factor  levels.  Finally,  it  was  noticed 
that  when  ilse  supply  temperature  as  measured  by  five  thermocouples  distributed 
over  the  center  of  the  supply  header,  was  above  room  temperature,  the  recovery 
temperatures  were  lower  than  expected,  and  vice  versa.  Surveys  of  the  flow  in 
the  tunnel  contraction  section  showed  strong  separation  so  that  it  was  likely 
that  the  air  in  the  center  of  the  header  was  not  the  air  that  passed  over  the 
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model.  The  correlation  wioh  supply  temperature  could  be  explained  if  it  was 
assumed  that  the  air  temperature  was  some  value  between  that  in  th»  center  of 
the  header  and  the  room  temperature.  Until  the  contraction  section  was  re- 
built, the  tests  were  ram  with  the  supply  temperature  set  within  i”  of  the 
room  temperature . 

TRANSITION  FIXED  ON  THE  CONF. 

Several  methods  were  tried  for  tripping  the  boundary  layer  on  the  cons, 
such  as  lampblack  and  sand.  However,  for  ease  in  i-unning  alternately,  in  the 
smooth  and  tripped  condition,  easily  removable  single  wire  trips  were  used. 

In  those  tests  the  wire  size  necessary  for  causing  turbulence  in  the  vicinity 
of  the  wire  was  larger  tnan  has  been  found  for  low-speed  measurements . For 
instance,  at  a Mach  Number  of  3.53,  at  an  equivalent  flat  plate  Reynolds  Number 
of  ?.15>  x 10^,  a wire  five  times  as  large  as  the  boundary  layer  displacement 

thickness  at  the  wire  did  not  cause  transition  anywhere  on  the  cone  surface. 

2 

According  to  Dryden,  at  low  speeds  at  the  acme  Reynolds  Number  a single  rough- 
ness element  .8  of  the  boundary  levvr  displacement  thickness  would  be  sufficient 
he  cause  transition  at  the  element. 

The  use  of  wire  trips  disclosed  30me  very  interesting  variations  of  the 
boundary  layer  recovery  temperature  as  influenced  by  the  size  of  the  wire  trip 
used.  While  these  variations  are  of  interest  themselves  and  should  be  the 
subject  for  another  investigation,  they  are  undesirable  for  this  investigation 
and  it  wa3  necessary  to  use  wire  sizes  for  which  they  were  not  present.  Figure 
12  presents  surface  temperature  data  at  the  same  Much  Number  and  Reynolds 
Number,  but  for  different  wire  sizes.  Shadowgraphs  for  these  three  cases  are 
shown  in  Figure  13.  The  smallest  wire  size,  .023  cm.,  produces  very  limited 
regions  of  separation  ahead  of  and  behind  the  wire  and  transition  occurs  about 
halfway  between  the  wire  and  the  model  shoulder.  For  the  median  wire  size, 
transition  starts  a short  distance  downstream  of  the  wire  and  again  the  separ- 
ation regions  ahead  of  and  behind  the  wire  are  relatively  limited.  For  the 
largest  wire  size,  separation  of  the  boundary  layer  occurs  well  ahead  of  the 
wire  and  transition  occurs  at  the  wire . There  is  a sharp  rise  in  the  surface 
temperature  ahead  of  the  wire,  and  a sharp  decrease  after  the  wire.  If  the 
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boundary  of  the  hoc’?!?.  fclleved  ih°  of  the  separated  region,  the  surface 

temperatures  would  be  expected  to  vary  as  they  do  because  of  the  changes  in 
surface  Mach  Number r However,  the  ienjper-aturss  in  the  folly' -devel  oped  tur- 
bulent region  on  the  cone  are  low,  and  they  remain  low  all  the  way  b^.ck  on 
the  cylinder,  only  gradually  approaching  the  torperature  data  with  the  smaller 
ware  triuo.  Static  pressure  measurements  along  the  model  surface  for  tvo  wire 
trips,  both  of  which  caused  ti-arsitlon,  are  cospared  with  the  smooth  model 
static  pressure  distribution  in  Figure  iiu  The  temperature  data  indicated 
that  t-b»  ,0«*  ea.  dia.  wire  was  too  large  to  be  satisfactory,  (see  Figure  12). 
The  pressure  changes  in  the  separated  region  ahead  of  and  dessstressi  of  the 
wire  trips  are  consistent  with  the  flow  photographs.  Cl  more  interest  is  the 
fact  that  the  surface  pressures  near  the  base  of  the  cone  are  only  ;»  •’lightly 
altered  by  the  wire  trips.  However,  the  pressure  model  was  also  fitted  with 
two  sets  of  three  total  head  tubes,  one  positioned  near  the  base  of  the  cone, 
and  one  placed  well  downstream  on  tb?  cylinder*  (C»s  Figure  U)  • The  total 
pressure  tubes  near  the  cone  surface  snowed  a loss  in  total  pressure,  presum- 
ably a result  of  the  shock  wave  system  associated  with  the  flow  about  the  wire 
trip.  For  the  ,023>  cm.  trip,  the  Mach  Number  decreased  by  about  ,Q6  M at-  the 
cone  base,  so  that  using  the  smooth  cone  surface  Mach  Number  for  calculating 
the  recovery  factor  3hcruld  have  resulted  in  the  c cupuycd  ? cs  i ng  i/w  by 

*00$ . The  Mach  Number'  decrease  at  the  cone  base  is  the  largest  with  the  .062 
cm.  wire  trip,  and  as  Figure  12  sherds,  the  computed  recovery  factors  are  lowered 
rather  than  raised.  The  greatest  Mach  Number  changes  on  the  cylinder  pr  educed 
by  the  wire  trips,  as  evaluated  frem  the  static  pressure  measurements  along  the 
cylinder  and  the  total  pressure  measurements  near  the  cylinder  sxirface,  would 
have  a negligible  influence  on  the  recovery  factor  c crrputations . Evidently,  the 
Mach  Number  changes  introduced  by  the  large  wire  tries  cannot  account  for  the 
effect  of  wire  size  on  the  temperature  data.  Possibly,  there  is  some  redistri- 
bution or  extraction  of  energy  involved  in  the  flow  changes  in  the  vicinity  of 
the  wire  vhiSh  persists  for  a very  large  number  of  boundary  layer  thicknesses. 
With  a two-wire  trip,  the  effect  is  even  more  pronounced.  On  the  basis  of 
these  experiments,  wire  sizes  were  established,  which  as  far  as  surface  temper- 
atures were  concerned,  produced  equilibrium  fully-developed  turbulent  boundary 
layers  just  ahead  of  the  shoulder. 
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A comparison  of  the  mcQ.crv  "s^oris  "or  the  smooth  models  and  or=  ease 
•where,  the  boundary  layer  was  turbulent  ahead  of  the  shoulder  is  shown  in  Fig- 
ure 11.  Shadowgraphs  of  the  flow  along  the  cylinder,  for  the  latter  case, 
show  clearly  that  the  boundary  layer  is  turbulent  from  about  8 cm . from  the 
model  shoulder  to  the  model  base.  See  Figure  15 . Also,  the  recovery  factors 
in  the  same  region  are  close  to  the  normal  levels-.  On  the  other  hand,  for 
about  5 S cm.  after  the  shoulder,  the  recovery  factors  appear  to  be  at  the 
laminar  levels.  Figure  16,  which  contains  data  for  tripped  boundary  layers 
obtained  over  a range  of  Reynolds  Numbers,  shows  that  the  data  at  various 
pressure  levels  lava  the  same  features.  Similar  measurements  at  K - 3*65  are 
shown  in  Figure  17.  At  Pn  - 100  cm.,  the  recovery  factors  appear  to  conform 
to  laminar  values  for  about  7.5  cm.  from  the  shoulder.  The  beginning  of  the 
temperature  "transition"  region  on  the  cylinder  is  difficult  to  pinpoint  be- 
cause of  the  relative  scarcity  of  the  temperature  measuring  positions.  At  any 
rate;  it  would  appear  to  be  more  than  2.5  cm.  from  the  3houlder  in  Figures  16 
and  17. 

Photographic  evidence  proved  to  be  difficult  to  obtain Figure  18  rep- 
resents the  best  shadowgraph  obtained  in  the  wind  tunnel.  The  foil. owing 
observations  can  be  made . The  turbulence  in  the  boundary  layer  ahead  of  th'. 
shoulder-  appears  to  disappear  in  the  expansion  at  the  shoulder.  A diffuse 
white  line,  sometimes  wavy  in  appearance,  can  be  seen  a short  distance  from 
the  shoulder,  and  then  turbulence  in  the  boundary  layer  appears  over  a range 
of  positions  which  agree  with  the  region  labeled  "transition"  based  on  the 
temperature  measurements . Other  wavy  laminar  boundary  layers , ju3t  prior  to 
transition  to  turbulence,  have  been  observed  under  conditions  where  the  laminar 
boundary  layer  was  much  thicker,  as  in  Figure  19. 

One  difficulty  with  getting  good  quality  wind  tunnel  shadowgraphs  is  the 

fact  that  the  density  in  the  test,  section  is  low.  A "ring"  was  machined  on 

the  nose  of  a 20mm  cone -cylinder  model  for  which  the  free  flight  Reynolds 

Number  correspoiided  to  one  of  the  wind  tunnel  Reynolds  Numbers,  and  t.he  model 

# _ 

was  fired  in  the  Ballistic  Research  T,aborat-r:ies  Aerodynamic  Range.  figure  20 
# The  firing  was  carried  out  under  the  supervision  of  L.  C.  McAllister* 
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is  a shadowgraph  of  the  free  flight  model.  Because  of  the  model  motion,  the 
turbulent  region  ahead  of  the  shoulder  is  more  difficult  to  distinguish;  how- 
ever, the  ring  sirs  chosen,  denied  arprwriately,  was  more  than  3 ufficient  to 
cause  transition  at  the  ring  in  the  wind  tunnel.  Just  after  the  shoulder  tnere 
is  the  characteristic  white  line  of  the  laminar  boundary'  layer  before  the  ap- 
pearance of  turbulence.  There  is  another  feature  in  this  photograph  which 
sugge>i»^  that  the  boundary  layer  is  leminar  a xt  ha3  been  observed  in  many  free 
flight  shadowgraphs  that  even  for  highly  polished  models,  (the  cylinder  of  the 
20mm  model  was  polished),  where  the  boundary  layer  is  turbulent,  there  is  a 
disturbed  pattern  in  the  surrounding  flow  field,  limited  cs  the  upstream  end 
by  the  Mach  wave  from  t<iie  beginning  of  the  turbulent  boundary  layer.  This 
disturbed  pattern  has  beer  assumed  to  be  a sound  field  originsting  from  the 
turbulent  boundary  layer.  In  Figure  20,  such  a pattern  is  visible  behind  the 
ring  on  the  cone,  (Mach  wave 3 are  also  visible;  they  arc  present  because  the 
model  surface  between  the  ring  and  the  shoulder  va3  not  polished),  and  in 
the  flow  field  around  the  cylinder.  The  Mach  wave  defining  the  leading  edge 
of  this  region,  about  the  cylinder,  extended  to  the  model  surface,  clearly 
intersects  the  model  downstream  of  the  shoulder  in  the  region  identified  as 
"transition*  on  the  cylinder# 


» 

3 i 


ANALYSIS 

These  experimental  results  can  be  explained  in  the  following  way. 

The  turbulent  boundary  layer  is  known  to  have  a thin  region  near  the 
surface,  which  is  identified  as  the  laminar  sub-layer,  where  the  direct 
viscous  terms  are  predominant.  Measurements  of  the  turbulent  velocity  fluctu- 
ation level  through  the  turbulent  boundary  layer  show  that  a peak  in  the  tur- 
bulence level  is  reached  at  the  edge  of  the  laminar  sub-layer.  This  high 


disturbance  level  apparently  prevents  the  laminar  sub-layer  from  growing,  so 
that  the  Reynolds  Number  based  on  the  thickness  of  the  laminar  sub-layer  is 
invariant. 


The  shaip  pressure  drop  at  the  model  shoulder  produces  a large  increase 
in  velocity  of  the  boundary  layer  fluid . wven  near  the  surface,  the  air  is 
accelerated  to  supersonic  speeds.  The  effect  of  the  velocity  increase  on 
the  turbulence  level  throughout  the  boundary  layer  ia  similar  to  t'ue  reduction 
in  turbulence  level  accomplished  by  the  contraction  section  of  a wind  tunnel. 
At  the  beginning  of  the  cylinder,  the  air  that  comprised  the  turbulent  bound- 
ary layer  Ahead  of  the  shoulder  is  the  non-uniform,  fairly  turbulent  "free 
stream"  for  a new  laminar  boundary  er  growing  from  the  beginning  of  the 
cylinder.  The  reduction  in  the  turbulence  level  of  the  boundary  layer  fluid 
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from  the  model  surface . Compared  to  the  usual  wind  tunnel  flow,  the  turbulence 


level  of  this  new  "free  stream"  is  high,  and  so  a transition  of  the  r.cv  laminar 
boundary  layer  to  a turbulent  boundary  la^’er  would  be  expected  to  occur  at 
relatively  low  Reynolds  Numbers  based  on  distances  along  the  cylinder  from  the 
shoulder.  The  only  turbulent  fluid  which  becomes  laminar  is  the  fluid  in  the 
new  laminar  boundary  layer,  "tlic  s &nc  process  that  occurs  in  tha  laminar  bound- 
ar-  layer  on  a model  in  the  usual  wind  tunnel,  since  the  wind  tunnel  air  is 
turbulent • 


i,  similar  acceleration  takes  place  when  there  is  a laminar  boundary  layer 
ahead  of  the  comer-  The  laminar  shear  flow  at  the  beginning  of  the  cylinder 
may  be  considered  to  be  a new  "free  stream"  for  a new  laminar  boundary  layer 
growing  from  the  shoulder.  The  evaluation  of  the  temperature  recovery  factors 
along  the  cylinder  must  take  into  account  the  fact  that  the  boundary  layer 
velocity  and  temperature  distributions  differ  importantly  from  the  velocity 
and  temperature  distributions  for  a boundary  layer  on  a constant  pressure 
surface*  (The  gradual  pressure  rise  along  the  cylinder  will  be  neglected  in 
the  computations.) 

In  the  following  sections,  estimates  will  be  made  for  the  boundary  layer 
conditions  on  the  cylinder  for  both  a laminar  and  a turbulent  boundary  layer 
at  the  base  of  the  cone.  For  the  case  where  the  boundary  layer  an  the  cone 
base  is  turbulent,  the  properties  of  the  turbid. ence  in  the  new  "free  stream" 
on  the  cylinder  will  be  examined  to  see  what  can  be  established  about  the 
probable  state,  laminar  or  turbulent,  of  the  new  boundaiy  layer  starting  from 
the  beginning  of  the  cylinder. 
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At  a number  of  points  in  the  analysis,  it  will  be  necessary  to  use  low- 
speed  data,  because  of  the  limited  nature  uf  supersonic  information- 


(a)  Turbulence  field  In  supersonic  boundaiy  layer. 

-5 

Recently,  Kovasznay  has  made  exploratory  velocity  and  t emperature 
fluctuation  level  measurements  in  a supersonic  boundary  layer.  Us  found,  in 
addition  to  the  velocity  fluctuation  field,  that  adjacent  fluid  masses  in  the 
boundary  layer  have  unequal  stagnation  temperatures;  i.e*  a "turbulent"  tempera- 
ture field  in  sn  Euler ian  system*  The  mora  voluminous  low-speed  data  are 
used  in  the  analysis,  but  a comparison  of  Kovaa^nay Ss  fluctuation  levels  and 
the  fluctuation  levels  obtained  frem  the  low-speed  measurements  will  be  made. 


\ 


3 


(b)  Mean  and  turbulent  properties  of  two-dimensional  supersonic  wakes. 


! 


The  author  was  unable  to  find  any  suitable  supersonic  data,  so  low- 
speed  data  have  been  used  exclusively, 

(c)  Influence  of  the  free  stream  turbulence  on  laminar  turbulent  tran- 
sition at  supersonic  speeds- 


The  available  information  in  the  literature  on  the  role  of  the 
turbulence  of  the  free  stream  on  the  state  of  the  boundary  layer  is  for  free 
streams  with  isotropic  turbulence.  When  the  boundary  layer  is  turbulent  at 
the  base  of  the  cone,  the  "free  stream"  at  the  beginning  of  the  cylinder  will 
be  anisotropic.  To  represent  the  "free  stream"  conditions  from  tbs  beginning 
of  the  cylinder,  this  anisotropic  field  Trill  be  represented  by  an  equivalent 

_ i 

isotropic  field  by  using  an  average  fluctuation  level  u and  a turbulence 
microscale  based  on  isotropic  relations . Along  the  cylinder,  turbulence  in 
the  "free  stream"  should  approach  isotropy.  As  will  be  seen,  the  interpre- 
tation of  the  results  is  not  sensitive  to  the  engineering  approximations 
that  are  made. 
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Lsting  solution  for-  the  compressible  laminar  boundary  layer  will  bo 
vimatc  the  effect  of  free  stream  turbulence  on  the  laminar  tvrWi!  ent 


transition  at  sunersonic  speeds. 


To  avoid  confusion  in  the  calculations,  the  following  notation  will  be 
adopted.  The  free  stream  Mach  number  distribution  along  the  cone-cylinder 
corresponding  to  the  static  pressure  distribution  will  be  designated  as  Mg. 
The  recovery  factor  r is  the  recovery  factor  computed  using  the  measured 
surface  temperature  Tg,  the  Mach  number  >L,  and  the  wind  tunnel  stagnation 
temperature  T n,  as  follows. 
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As  a new  boundary  layer  grows  along  the  cylinder,  its  edge  at  ary  axial 
position  is  at  some  Mach  Number  less  than  Mg . At  ary  axial  position  the 
value  of  this  ‘’effective"  free  stream  Mach  Number  will  depend  on  the  new 
boundary  layer  thickness  at  that  position  and  tne  properties  of  the  n free 
stream"  flow.  Primes  will  be  used  to  identify  the  conditions  at  the  edge 

I 

of  the  "new*1  boundary  layers.  Hence,  r would  be  obtained  as  follows. 

. - .2- 

T — T • 1 4-  Of r,  i 

i _ "o  'SI1  *“T5  I 


i _ *o  SI  *“T5  I 
r - x r— , — o 

l0  L_  ...  • J 

•2hs 

Tc  - stagnation  temperature  of  flow  at  edge  of  "new" 
boundary  layer 

Mg  = Mach  Number  at  edge  of  "new"  boundary  layer 
Tg  « Surface  temperature 

rT  - Laminar  recovery  factor  = .35 
x» 

rT  ■ Turbulent  recovery  factor  =»  .88 

t i 

VThila  Tq  and  Mg  are  insufficient  to  establish  properly  the  temper ature 
recovery  factors  along  the  cylinder,  the  above  expression  is  considered  to 
be  a useful  approximation. 
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As  shewn  in  Figures  16  and  1?,  the  bonporature  data  at  different 


pressure  levels  all  had  the  same  general  characteristics  . The  numerical 


calc- ula tiens  that  arc  made  in  the  following  sections  are  for  a supply 
pressure  * 180  cm.  at  M = 3o02,  which  is  considered  to  be  a representa- 


tive test  condition.  Data  for  a turbulent  boundary  layer  at  the  cone  base 
are  for  a foUl  cm.  boundary  layer  trip. 


TURBULENT  BOUNDARY  IAXiiK  .V-  akSiL  OF  CONE 


The  shadowgraphs  indicate  that  the  boundary  layer  with  the  .OlfL  cm. 


trip  is  turbulent  for  about  UO  boundary  layer  thickne33es  before  it  retches 

— A/Snn  W.  ma  A . « n nj  M — 1—  n -I  rt f A'  rtVI  ^ A'  r~-~  ■!  A tl  O A Vt*.  /A  A C~ 


the  cone  bo.se-  At  law  speeds,  Klebanoff  ard  DieU/*  found  that  badly  dis- 
turbed turbulent  boundary  layers  had  re-ostablished  the  equilibrium  or 
constant,  prssc'xe  surface  velocity  profiles  and  t-irbulont  fields  after 
traveling  5b  30  boundary  layer  thicknesses  - Further  evidence  that  equilibrium 
conditions  are  closely  approached,  if  not  reached,  at  the  cone  base  is  the 
fact  that  the  tempe rature  recovery  measurewaats  near  the  cone  base  correspond 
to  those  shown  in  Fig-ore  7.  Available  measurements  Indicate  that  the  velocity 
profiles  in  supersonic  boundary  layers  are  similar  to  those  in  low-speed  bound- 
ary layers,  provided  the  Reynolds  Numbers  are  comparable . It  is  therefore 
assumed  that  the  boundary  layer  velocity  can  be  represented  as: 


, „ ,1/7 
" . It) 


Since  transition  is  produced  by  a roughness  element,  it  is  net  possible 
to  compute  the  boundary  layer  thickness  at  the  cone  base.  The  shadowgraphs 
show,  as  would  be  expected,  different  apparent  boundary  layer  thicknesses  for 
different  3ize  roughness  elements,  at  different  pressure  levels,  and  different 
Mach  Numbers. 


A large  fraction  of  the  turbulent  boundary  layer  has  an  intermittent 
character,  which  is  three-dimensional  in  nature.  Low-speed  measurements'’’ ^ 


shew  that  there  is  ?.  sharp  boundary  between  the  turbulent  fluid  in  the  bound- 
ary layer  near  ths  wall  and  the  main  flow  0 This  boundary  is  irregular  and  in 
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a typical  case  may  extend  from  £ = .U  to  ? 3 1.2*  with  h mean  position  at 

^ “ ,?8,  where  6 is  the  conventionally  defined  boundary  layer  thickness ; 

In  one  case,  a wave  length  from  oeak  to  peak  of  the  turbulent  boundary  of 
26  was  observed  using  hot-wire  equipment. 

The  shadowgraphs  at  high  speed  of  turbulent  boundary  layers  and  wakes 
show  the  3anvs  sharp  boundary  between  the  turbulent  and  apparently  non* 
turbulent  fluid.  The  strong  optical  effect  is  presumed  to  be  due  to  the 
temperature  spottiness  of  the  turbulent  fluid.  Since  the  path  length  of  tbs 
light  rays  in  the  boundary  layer  is  usually  several  boundary  layer  thicknesses, 
the  shadowgraphs  do  not  show  the  large  variations  in  the  position  of  the  edge 
of  the  turbulent  boundary  layer  that  are  found  in  the  low-speed  measurements. 
However , the  boundary  has  a sharp'l2r  irregular  appearance.  If  the  path  length 
of  the  light  through  the  boundary  layer  was  very  large  compared  to  6,  then  it 
could  be  argued  that  the  boundary  in  the  shadowgraph  should  represent  1 .25 . 

This  requirement  is  not.  fulfilled  for  the  test  model,  and  since  peak  excursions 
of  the  turbulent  fluid  may  not  be  at  the  top  cone  element,  it  is  difficult  to 
establish,  with  asy  precision,  the  boundary  layer  thicknesses  from  the  shadow- 
graphs . On  the  other  hand,  it  seems  reasonable  to  believe  that  the  photographic 
boundary  lies  somewhere  between  .786  and  1.26.  The  simplest  assumption  is  to 
set  the  apparent  thickness  on  the  photographs  equal  to  6, 

An  appropriate  picture  of  the  boundary  layer  velocity  and  temperature 
profiles  can  be  obtained  by  relating  the  local  friction  coefficient  to  the 
boundary  layer  thickness , For  a compressible  boundary  layer  with  zero  axial 

2 /40  Jq  Q 

pressure  gradient,  we  have  tr^-  “ pgUg  or  -y"  = 3x  w*lere  the 

local  friction  coefficient. 

Now,  for  incompressible  fluids,  in  the  range  of  Reynolds  Numbers  where 
the  velocity  profile  can  be  represented  by  a i/?ih  power  law,  ws  can  express 
the  local  friction  coefficient  in  terms  of  the  length  Reynolds  Number  as 


£*• 

4 


*r««' " . 


Rocently,  a number  of  experiments  have  been  performed  it  supersonic 
speeds,  in  which  the  local  friction  coefficient  has  been  compared  with  the 
incompressible  local  friction  coefficient  at  the  same  Reynolds  Number.  The 
compress ible  local  friction  coefficient  can  be  represented  as  follows: 

0 = -°^T/5  f(M) 

« ^x  ' 


where  f (M)  is  given,  for  instance,  ii 
Then  £ - K f (M)  or  |-  f (M) 


- ’TA  r A <!•«!  90 

x,  »»>v  x..  r - — ✓ ■ 


Re  ~ 
x 


f)  -5. 

To  compare  ~ compressible  and  — Incompressible  at  the  same-  x Reynolds 

0 X 

Number,  the  ration  g is  required.  Assuming  the  Prandtl  Number  equals  1,  and 
no  heat  transfer  at  the  surface,  it  can  be  shown  that 


JL  _ 

e f F(n)  Li  - F(n)  1 

* o 1 + ^5“  C1  * F (yP] 
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Convenient  tabulations  of  the  value  of  this  ratio  for  various  Mach  Numbers 

and  various  power  law  representations  of  the  velocity  distribution  are  given  in 

NASA  TN  2337.  If  f » 9(M) 

then  6C  - °lJ»  O(M-O)  f(M) 
then  ^ 

The  factor is  shown  in  Figure  21,  and  it  appears  that  for  a given 
free  stream  Reynolds  Number,  the  ration  ~ varies  less  than  with  Mach 
Number,  Approximately,  then  6C  ■ 6^c,  and 

•Oli$  f(M) 

c*.  it-™ 

°c 

Meajrurement  of  the  shadow  graph  gives  6 - ,06  cm.  at  the  ba36  of  tha  cone, 
so  that  Seg  - 1.?  x lO*4  Cf(H.  - 1.69)  - 2.9  x 10”3 


Using  the  measured  wall  temperatures 

<s 

30  that  UT. « 2oh3  x 10-5  cm/sec 


XW  , 

. m-  * ls»\i 

As 


i >•  v 

-„nsrc  V r - - ' — the  friction  velocity, 

' Py  / 

Finally,  at  the  wall  ( ) » 1,20  sec”^ 

' / W 

The  extent  of  the  laminar  sub-layer  is  a matter  of  definition . Some 
UTy 

authors  choose  = 10  as  the  edge  of  the  layer  since  that  is  where  the 
direct  viscous  stress  and  the  turbulent  shear  stress  are  approximately  equal, 

H_y 

Others  choose  - ^ ■ 30,  where  the  direct  viscous  stress  is  negligible  and  t-he 

mean  velocity  profile  starts  to  follow  the  logarithmic  law.  From  the  wall  to 
U~y  ^ tfTy 

— — = 10  the  velocity  profile  is  linear  and  the  region  = 10  to  — * 30 

represents  a transition  region  to  a logarithmic  mean  velocity  profile.  The 
construction  of  the  desired  velocity  profile  was  simplified  as  follows: 

U = ky  near  the  wall  where  (•!—)  =*  1.20  sec  * 77  = { *’or  J’ost 

\dy/w  us  n°/ 

of  the  boundary  layer. 


The  linear  variation  of  U with  y was  continued  until  the  profile  based 
on  the  l/7th  power  law  was  intersected  so  that  no  attempt  was  made  to  introduce 
the  correct  distribution  for  the  intermediate  region  where  the  laminar  shearing 
stresses  and  the  turbulent  shearing  stresses  are  Doth  important . Figure  22 
show3  the  Mach  Number  distribution  assuming  Pr  si.  (Because  of  this  assump- 
tion, the  Mach  Numbers  near  the  wall  are  a few  pnr  cent  in  error . ) 


The  fact  that,  the  acceleration  of  the  flw:  around  the  corner  occurs  in  a 
relatively  short  distance  makes  possible  a simple  computation  of  the  boundary 
1 aver  crofile  after  the  expansion.  The  mcxientum  equation  for  the  x direction 
of  the  low-speed  turbulent  boundary  layer  may  be  written  as 

p ( u H + v f ) " h [u  % + ~ sr 

where  u ^ is  the  laminar  shearing  stress  and  puv  is  the  Reynolds  shear  stress, 
cy 

and  the  forces  are  expressed  as  forces  per  unit  mass  of  tho  fluid.  The  shear 
stress  is  a maximu...  at  the  wall,  and  decreases  steadily  through  the  boundary 
layer  to  close  to  aero  at  y ■ 6 = 


11 


Then,  approximately,  the  viscous  i'orco  per  unit  mas  of  the  boumary 

f ___ 

layer  fluid  is  ^ « In  a typical  '"rrbulsrrt,  boundary  layer ^ - -3  fr«s 


the  wall  to  £ » .8,  failing  off  rapidly  to  z^ro  between  ^ - „8  and  ^ =•  1.2 
(the  maximum  extent  of  the  turbulent  spikes) , The  effect  of  the  acceleration 
at  the  corner  on  the  correlation  coefficient  ~r$T  i3  not  known*  It  will 
develop  in  a later  section  that  the  product  u'v*  after  the  comer  is  reduced 
to  about  -i|2  of  it3  value  ahead  of  the  comer,  so  that  even  if  -22—  after  the 

’2'T* 

corner  was  increased  to  the  maximum  possible  value  of  one , the  average  uv 
would  be  reduced  in  the  pressure  drop  at  the  shoulder.  Since  the  turbulent 


boundary  layer  thickness  increases  almost  by  a factor  of  3,  (See  Page  28)  it 

rw 

appears  that  -^=  throughout  the 
just  before  the  acceleration. 


appears  that  rg  throughout  ths  turbulent  boundary  layer  will  be  a maximum 
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Referring  to  the  x momentum  equation  for  the  boundary  layer,  for  a very 

- 3PS 

rapid  and  substantial  acceleration  the  term  i . which  is  the  pressure 

force  per  unit  mass  of  the  boundary  layer  fluid,  may  be  very  large  compared 

v W ^ W 1 ~PS 

to  . If  can  be  neglected  compared  witn  — , then  the  Mach  number 

profile  after  the  corner  can  be  obtained  in  a simple  way.  The  relative  in- 
fluence of  the  viscous  and  pressure  terms  can  be  evaluated  by  considering 
their  contribution  to  the  change  in  kinetic  energy  pet  unit  mass  of  the  boundary 
layer  fluid  in  the  acceleration  at  the  shoulder » The  acceleration  length  is 
greatest  for  the  outer  streamlines,  where  it  may  be  from  ? to  10  6,  but  small 

for  the  inner  regions,  where  it  should  be  of  the  order' of  5V  Taking  the  acceler- 

Ty  ' U 

ation  length  as  6,  AK.E^.  due  to  Viscous  Forces  ~ ^ x 6»  — 

* n 

A K.E.p  due  to  Pressure  Forces » ^XS 

P 

For  -'j i increase  in  Mach  number  from  1.89  to  3.13 
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viscous  forces  can  be  neglected. 


The  Mach  Number  for  each  streantube  after  the  expansion  can  then  be 

easily  computed  from  the  Mach  Number  of  the  streamtube  ahead  of  the  corner, 

assuming  that  the  flow  in  each  streamtube  goes  isentropically  from  Ps  to  P_  . 

‘l  '2 

The  dependence  of  the  Mach  Number  on  the  distance  from  the  wall  after  the 
acceleration  i3  obtained  by  a numerical  step  by  step  calculation  for  the  bound- 
ary layer  mass  flow  from  the  wall . Figure  22  shows  the  results  of  these  com- 
putations. The  air  that  was  almost  stationary  ahead  of  the  shoulder,  adjacent 
to  the  wall,  is  accelerated  to  almost  M ■ 2 after  the  shoulder.  Of  course  the 
conditions  right  in  the  vicinity  of  the  wall  are  not  properly  represented  by 


the  approximation,  since  the  viscous  forces  there  cannot  be  neglected.  At  the 
wall,  the  momentum  equation  reduces  to  ^ ^jand  for  the  real  flew 
the  Mach  Number  will  drop  rapidly  to  zero  in  a narrow t region  near  the  wall. 


The  expansion  of  the  streamtubes  results  in  an  increase  of  the  boundary 
layer  thickness  6 from  6 


•08  cm . to  5 * 


.23  cm. 


The  boundary  layer  fluid  at  the  beginning  of  the  cylinder  can  be  viewed 
as  a "new”  free  stream  for  the  development  of  a new  boundary  layer  downstream 
from  the  shoulder.  Its  Mach  Number  ranges  from  M ■ 2 at  the  wall,  tc  M = 3»13 
at  6 - .23  cm.,  beyor-d  which  distance  the  free  stream  is  approximately  uniform. 
This  "new”  free  stream  should  not  retain  the  Mach  Number  distribution  which  has 
been  calculated  for  the  start  of  the  cylinder.  Turbulent  mixing  should  alter 
the  Mach  Number  distribution  and  should  cause  a spread  of  the  mixing  region 
into  the  uniform  outer  flow.  If  the  shear  along  the  cylinder  surface  were  zero, 
then  the  turbulent  mixing  region  would  correspond  to  half  of  a two-dimensional 
turbulent  wake.  Presumably,  the  Mach  Number  distribution  which  has  been  com- 
puted for  the  flow  after  the  comer  would  approach  the  equilibrium  distribution 
for  a turbulent  wake.  In  the  real  flow,  the  existence  of  shear  at  the  wall 
requires  the  growth  of  a new  shear  la>er  starting  at  the  shoulder  which  could 
be  either  a new  laminar  boundary  ?ayer  or  a new  turbulent  boundary  layer.  Tne 


2d 


growth  and  properties,  such  as  temperature  distribution,  of  either  of  these 
new  boundary  icy  ore  ie  complicated  by  (a)  the  variable  Mach  Ifumbcr  distri- 
bution of  the  free  stream,  and  (b)  the  variable  stagnation  enthalpy  or  tne 
free  stream. 


r'here  is  a redistribution  of  the  stagnation  enthalpy  in  the  boundary- 
layer  fluid  ahead  of  the  corner-,  with  a decrease  near  the  wall  and  a com- 
pensating increase  in  the  outer  region  of  the  bo-uidary  layer.  The  miring 
after  the  turn  may  cause  a further  change  in  the  enthalpy  distribution. 

The  new  boundary  layer  growth  will  be  estimated  in  two  different  "free 
streams" . 

1.  The  boundary  layer  Mach  Number  distribution  as  calculated 
for  the  beginning  of  the  cylinder.  (See  Figure  22.)  This  "free  stream" 
will  remain  unchanged  along  the  cylinder,. 

2.  A two-dimensici-ial  wake  with  the  model  surface  corresponding 
to  the  center.  The  wake  width  will  increase  and  the  speed  defect-  near  the 
surface  will  decrease  with  Increasing  distance  along  the  cylinder. 

EQUIVALENT  two-dimensional  wake 


First  it  will  be  assumed  that  friction  at  the  model  surface  is  zero, 
so  that  the  model  surface  corresponds  to  the  wake  center.  The  influence 


of  friction  at  the  wall  can  be  japarately  evaluated.  The  asymptotic  rep- 

7 

resentetion  for  the  mean  velocity'  distribution  in  a two-dimensional  wake 

|—  3/2 — I 2 

is  tt  - 1 -( %—')  where  U is  the  maximum  v elocity  deflect 

Vz  L W J max 

oc curing  at  the  wake  center  and  y^  equals  half  of  the  wake  width  as  defined 
by  the  above  equation® 


The  two  conditions  necessary  to  obtain  and  yc  at  the  beginning  of 
tho  cylinder  are  chosen  so  that  the  wake  and  the  boundary  layer  are  equivalent. 
Since  the  momentum  of  a wake  does  not  vary  with  distance  downstream,  the 
momentum  derect  OX  X» njg  Whuri o equal  to  the  momentum  defect  of  the  bound- 

ary layer, 
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It  also  seems  reasonable  to  sst  the  total  mass  of  fluid  in  the  equivalent 
wake  equal  to  the  mass  of  fluid  in  the  boundary  layer* 
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The  equation  for  the  displacement  thickness  5 , for  PK  = 1 and  zero  heat 
transfer  is  ^ r ^ 

Using  these  two  equations,  and  numerically  evaluating  the  integrals  for 
U * 

and  6.^,  and  yQ  can  be  found  in  terms  of  00^  and  50  . 


For  the  equivalent  wake  at  the  begiimlng  of  the  cylinder, 
II 


max 

r— 

'0 


- .08U  and  7()  = .99  60J  > 


To  calculate  the  d ownstream  development  of  this  equivalent  wake  it  is 

7 

necessary  to  obtain  an  effective  starting  distance.  From  Goldstein, 


max 
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Hence,  an  effective  starting  distance  is  x - 3*5  cm.  along  Jin  equivalent 
cylinder • 


Q 

Townsend0  has  measured  the  turbulence  in  a low-speed  two-dimonnionel 
vake  and  has  found  that  the  turbulence  Reynolds  Number  RaT^  ^ RecyT  where 

P-  « , -s-  SA£  D --  » cylinder  diameter  and  He,™  * _ . ''LM2 

cyl  v eff  T v 
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For  a distance  £±  cw  who  cylinder  uf  souui;  100  e-y  Under  diameters  « n ~..3 

utnax 

independent  of  the  Reynolds  Number . Setting  Re„  * j P-e^^  results  in 
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that  tho  flow  would  not  closely  approach  the  wake  distribution  in  one  or  two 
cun.  Nance.  the  mean  flow  distribution  probably  lies  between  the  boundary 
layer  distribution  at  the  beginning  of  the  cylinder  and  the  wake  profile* 

The  effective  wake  width  would  increase,  proceeding  dewnstreem,  pro- 

■)  /J5 

portion  a!  to  x~'  , where  x - 3 *5  cm . plus  the  distance  from  the  shoulder « 


The  supersonic  laminar  boundary  layer  on  a flat  plate  grows  according 
to  the  following  law. 

3x  ” F Tf — ^ where  ^ - h(l^)  and  ■ g(Ji-) 

ih2%> 6 Vs 


As  an  engineering  approximation,  the  boundary  layer  growth,  was  numeri- 

* « ! 

cally  computed,  step  ty  step,  using  appropriate  values  of  g(KL  Y and  h(Hj  ) 
?_+.  each  step  of  the  process,  for  the  Mach  Number  distribution  of  the  bound  - 
ary layer  after  the  comer  and  for  the  equivalent  wake  Mach  Number  diatr-i- 
bution  along  the  cylinder. 

There  was  little  difference  between  the  results  for  the  boundary  layer 
and  wake  distributions,  so  that  apparently  & knowledge  of  the  exact  Mach 
Number  distribution  of  the  "free  stream"  is  not  necessary.  The  new  laminar 
boundary  layer  growth  in  the  equivalent  wake  is  shown  in  Figure  23. 

The  estimates  for  the  development  of  the  wake  "free  stream"  along  the 
cylinder  have  been  made  with  the  friction  at  the  surface  of  the  cylinder  set 
equal  to  zero.  As  a new  laminar  or  turbulent  boundary  layer  grows  from  the 
beginning  of  the  cylinder,  the  "free  stream*  fluid  is  consumed  so  that  the 
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momentum  dcfact  of  the  wake  ’’free  stream”  is  steadily  decreased  and  its 
development  should  be  altered.  The  following  calculations  made  for  conditions 
U.5  cm,  from  the  shoulder  irkiicate  that  an  far  as  the  calculation  for  the 
growth  of  the  new  laminar  boundary  layer  is  concerned,  there  is  little  error 
in  using  the  mean  velocity  distribution  for  the  wake  "free  stream"  obtained 
by  setting  the  surface  friction  equal  to  cero. 


At  x » U.5  cm.  the  wake  "free  stream”  width  yQ 
new  laminar  boundary  layer  thickness  6 ^ .2yn . 


x 


10  "3 


cm.  and  the 


At  x ■ k*5  cm,  a new  equivalent  wake  is  determined  by  considering  the 
edge  of  the  new  laminar  boundary  layer  to  be  toe  center  of  the  new  equivalent 
wake  and  by  equating  the  momentum  defect  and  mass  flow  for  the  new  equivalent 
wake  yQ  and  the  truncated  wake  from  y - 6 to  y - y,.,  (see  Figure  2U) . 


There  results  yQ  + 5«y0 

so  that  the  outer  edge  of  the  new  equivalent  wake  is  about  the  same  as  for 
the  truncated  wake, 

i 

U 

max  i . t 

Also,  = 4o  x iu  la  very  close  to 

S 

TT 

- e ^2 

~ - U.6  x 10”  , the  mean  velocity  in  the  truncated  wake  at  y - 6. 


GROWTH  OF  i!EW  TURBULENT  boundary  layer 


It  was  previously  chcvn  that  for  a given  Re^.,  that  ^ varied  very  little 
with  Mach  Number,  so  that  it  is  sufficient  to  take 

where  the  primes  denote  the  conditions  in  the 

"free  stream"  at  the  edge  of  the  new  turbulent  boundary  layer.  A3  for  the 
new  laminar  boundary  layer,  the  new  turbulent  boundary  layer  growth  is  not 
sensitive  to  small  changes  in  the  "free  stream”  flow.  Figure  23  shows  the 
results  for  the  growth  of  a new  turbulent  boundary  layer  from  the  beginning 
cf  the  cylinder  in  the  wake  "free  stream’* . 
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lakenah  boundary  layer  at  bask  of  cone 


i.  &pprcsd==tios  ussd.  for  obtaining  the  effect  of  the  pressure  drop 
on  the  turbulent  boundary  layer  mean  velocity  distribution  should  be  just 
as  valid  when  the  boundary  layer  at  the  base  of  the  cone  is  laminar.  At 
- 180  cm.,  '•'*  - 1.8?,  the  laminar  boundary  layer  thickness  at  the  base 
of  the  cone  is  2.5  x 10*”~  cm.  Figure  2p  snows  Hie  cCBiputed  velocity  pro- 
files before  and  after  the  corner.  (For  simplicity,  a linear  profile  was 
used  ahead  of  the  comer.)  In  this  case,  the  boundary  layer  thickness  is 
not  increased  by  as  large  a factor  as  for  the  turbulent  case . Since  a 
much  greater  extent  of  the  laminar  boundary  layer  ahead  of  the  corner  is 
subsonic,  the  not  expansion  of  the  streamtubes  is  less. 

However,  as  compared  with  the  turbulent  case,  there  is  an  important 
difference.  In  the  turbulent  wake,  the  turbulence  mechanism  provides  some 
interchange  across  the  whole  extent  of  the  wake  since  there  will  be  some 
eddies  of  a size  comparable  to  the  wake  width.  For  the  laminar  shear  flow, 
the  inhere  hang  a depends  on  the  molecular  transport  mechanism  and  has  a 
local  behavior.  Consider  the  following  cases.  A laminar  boundary  layer 
grows  from  the  leading  edge  of  a flat  plate  at  rest.  (See  Figure  3U » ) 


Between  x,  and  a movable  surface  is  inserted.  Now  the  boundary 
layer  growing  from  the  beginning  of  the  plate  represents  the  propagation  of 
significant  information,  by  viscous  forces,  of  the  presence  of  the  wall. 

An  element  of  fluid  at  A has  had  a slight  change  in  its  motion  because  of 
a change  in  the  flow  potential,  but  it  has  not  yet  been  significantly  af- 
fected by  viscous  forces , The  fact  that  the  surface  is  moving,  starting 
at  x^,  propagates  along  some  surface  such  as  the  one  shown.  This  infor- 
mation should  be  transmitted  away  fr  Oui  wiic  plate  more  rapidly  by  viscous 
forces  than  for  the  initial  growth  of  the  boundary  layer  from  the  plate 
leading  edge  because  of  th**  1 'wnr  velocities  near  the  wall.  The  fluid 
element  at  B has  not  learned  of  the  change  in  boundary  condition  at  x, 
and  decreases  its  velocity  at  the  same  rate  as  if  the  wall  condition  at 
x^  had  been  unchanged. 


With  the  velocity  distribution  after  the  corner  shown  in  Figure  25, 
the  fluid  elements  of  the  new  ires  stream  should  move  roughly  0.3  if  they 
belonged  to  another  laminar  boundary  layer  of  much  greater  thickness . 

The  profile  after  the  corner  is  almost  linear  and  might  be  imagined  to 
be  the  outer  reaches  of  a boundary  layer  with  a linear  profile  growing 
cn  a surface  well  below  the  r-eal  surface. 

Since  the  boundary  layer  growth  decreases  with  boundary  layer  thick- 
ness, the  new  laminar  boundary  layer  should  engulf  the  old  laminar*  layer 
about  2 o cm.  from  the  corner.  Probably,  some  further  travel  downstream 
would  be  required  to  reach  an  equilibrium  profile.  (Since  for  air  the 
temperature  boundary  layer  is  of  the  same  order  as  the  viscous  layer,  what 
has  been  said  for  sudden  changes  in  shear  at  the  surface  would  hold  equally 
well  for  sudden  changes  in  surface  temperature.) 

TEMPERATURE  RECOVERY  FACTORS  ALONG  THE  CYLINDER 


In  the  discussion  of  the  experimental  results  the  temperatures  on  the 
cylinder  were  described  in  relation  to  the  "usual"  laminar  or  "usual"  tur- 
bulent, levers#  it  is  crear  ircm  one  i urOg’’-* 9. v* that  the 
temperature  recovery  factors  for  the  new  laminar  or  turbulent  boundary 
layers  on  the  cylinder  would  differ  from  the  constant  pressure  surface 
levels • 

» ! 

Calculations  will  be  made  assuming  that  in  terms  of  m and  T0  , 

II  O w 

rT  =»  -85>  and  r™  = .88 . (The  1 imitations  of  this  assumption  were  pre- 
h 1 

viously  noted.)  The  recovery  factors,  r,  along  the  cylinder  for  the 
various  cases  can  then  be  obtained  in  terms  of  Mg  and  Tq. 

As  there  is  heat  transfer  between  the  air  and  the  model  for  about 
one  centimeter  ahead  of  and  downstream  of  the  model  shoulder,  the  measured 
surface  temperatures  in  the  corresponding  regions  are  not  recovery  temp- 
eratures. The  recovery  factor  calculations  for  the  boundary  layers  on 


the  cylinder  wilr  not  Lake  this  heat  transfer  Into  account  s«  that  the 
measured  tempera  Lures  in  the  vicinity  of  the  shoulder  should  not  agree 
with  the  calculations. 

(a)  Laminar  at  base  or  cone. 

For  rT  - .85  at  the  cone  base i,  ■ .937  where  Me  = 1.89,  so 

L T.  ~ 

u » 

that  for  the  fluid  adjacent  to  the  wall  ahead  of  the  shoulder  TQ  * T^.  . 

A new  laminar  bcundaiy  layer  at  the  beginning  of  the  cylinder  would  have 
KL  - 1.91  and  tJ  » .937  TQ  and  with  rT  ' =>  .85,  r * .82  based  on  JL  - 3.13 
and  a stagnation  temperature  T_. 


(b)  Turbulent  at  base  of  cone. 


For  r^,  * .88  at  the  cone  base,  1^  * .95  Tq  so 


that  for  the  fluid 


adjacent  to  the  wall  ahead  of  the  shoulder  = T^  . If  the  new  boundary 

» , « 

layer  after  the  corner  is  laminar,  then  with  rr  » .85,  Mg  « 1.91,  and 

Tq  - .95  T0, - r * .834  at  the  beginning  of  the  cylinder  based  on 

% - 3-13  • I T.. 


i 

i 


If,  as  previously  estimated,  the  new  laminar  boundary  layer  engulfed 
the  old  laminar  boundary  layer  by  2.5  cm.  f ran  the  corner,  then  from  2.5  cm. 
downstream  the  free  stream  for  the  growing  boundary  layer  would  be  that 
obtained  from  the  surface  pressure  variation.  Figure  11  3hws  that  the 
recovery  factor  r has  only  reached  a value  of  .84  by  2.5  cm.  from  the  corear, 
and  appears  to  level  out  at  a somewhat  higher  value  than  expected.  The 
reason  for  this  behavior  is  not  known;  it  may  be  that  after  engulfing  the 
old  boundary  layer,  the  new  laminar  layer  approaches  the  usual  equilibrium 
temperature  distribution  fairly  slowly.  Certainly,  it  does  not  seem  likely 
that  the  temperature  distribution  at  the  point  of  engulf ment  is  the  equi- 
librium distribution  since  the  old  boundary  layer  suffered  a pronounced 
energy  redistribution  ahead  of  the  comer,  and  this  should  effect  the  distri- 
bution of  energy  in  the  new  boundary  layer. 
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Kova3anayJ  found  that  the  total  enthalpy  was  constant  within  his 
measuring  accuracy,  (approximately  1%)  across  a supersonic  boundary  layer 
at  a riaeh  Number  of  1,7?,  at  least  between  £ - .05  and  2-1,  The  measure” 
meats  covered  the  supersonic  portion  of  the  boundary  layer.  Evidently,  the 
loss  in  enthalpy  occurs  vary  near  ths  wall.  ^ = ,1  in  the  turbulent  layer 
ahead  of  the  corner  corresponds  toy  » 1,6  x lO^*"  cm.  after  the  comer.  The 
new  laminar  layer  would  have  grown  to  this  thickness  within  about  *2  cm,  from 
the  shoulder-  Therefore,  except  for  conditions  ^iust  after  the  shoulder,  the 
"frea  stream"  enthalpy  would  be  approximately  the  stagnation  enthalpy  of  the 
tunnel,  but  the  "free  stream"  Mach  Number  would  be  reduced.  Because  of  the 
reduced  "free  stream"  Mach  Number,  the  apparent  recovery  factor,  r,  would  be 
greater  than  .85. 


1 f 

If  the  new  boundary  layer  is  turbulent,  then  rT  - .88,  * 1.91, 

Tq  “ .95  T0,  and  r ■ ,85  at  the  beginning  of  the  cylinder. 

The  apparent  recovery  factor,  r,  along  the  cylinder  for  the  turbulent 
boundary  layer  would  be  greater  than  r * .88.  Figure  27  compares  the  cal- 
culated apparent  recovery  factors  and  the  experimental  results-  A ~\$  dif- 
ference between  the  actual  and  assumed  stagnati  Oi>  te?iiip«i"fi-biit'HM  Of  th«  ”irKM 


stream",  due  either  to  a variation  ahead  of  the  comer  or  a redistribution 
after  the  corner,  would  correspond  to  a change  of  the  predicted  recovery 
factor  levels  for  the  new  laminar  and  turbulent  boundary  layers  of  about  .015  • 


The  course  of  the  recovery  temperatures,  as  shown  in  Figures  16  and  11, 
and  perhaps  more  clearly  at  M = 3-55  in  Figure  17,  is  a rise  from  the  shoulder 
flattening  out  at  2 - 7 cm.  from  the  shoulder  and  then  a fairly  steep  rise  to 
much  higher  values . 

The  conclusion  that  the  sharp  rise  in  the  experimental  temperatures  on 
tha  cylinder  is  due  to  a transition  from  a "ns w"  laminar  boundary  layer  to 
a "new41  turbulent  boundary  layer  is  believed  to  b9  consistent  with  these 
approximate  calculations . 
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TRANSITION  BY  STRONG  TURBULENCE 


At  low  speeds  the  3_ oi. wi i liuiu  ramma*.  v:-  "cLi. w-uiouv  xavv  Lh  kno=xj 

to  depend  on  the  turbulence  of  the  wind  tunnel  air3tream.  Taylor10  pro- 
posed a suitable  parameter  for  correlating  the  transition  results  when  the 
turbulence  of  the  stream  was  relatively  large » In  the  Karman-Pohlhausen 
method  for  calculating  the  behavior  of  laminar  boundary  layers  in  the 
presence  of  pressure  gradients,  the  boundary  layer  profile  at  any  position 
depends  on  the  local  boundary  layer  thickness  and  the  local  pressure  grad- 
ient, and  is  defined  in  terms  of  a parameter  _A.  where 


(It  is  known  that  for  rapid  changes  in  pressure  the  local  profile  depends 
on  the  boundary  layer  history  and  cannot  be  described  in  this  simple  way.) 
Taylor  s uggested  that  in  the  presence  of  strong  turbulence , the  instan- 
taneous pressure  gradient,  due  to  the  turbulence,  caused  separation  of  the 
laminar  boundary  layer  with  resulting  transition.  Taking  the  Pohlhausen 
treatment  as  a useful  approximation,  in  spite  of  its  known  limitations,  and 
estimating  the  root  mean  square  value  of  the  sp3.tial  pressure  gradient  in 
isotropic  turbulence,  that  is 


up  u’ 


.2  .2  Xn  ■ Lagrangian  microscale 

ap’  P*P  u1  op  U' 

„ •*  — ~ —nr X ■>  Longitudinal  microscale 

^ n x 

*•  u’  ■ root  mean  square 

fluctuation  velocity 

then  the  Pohlhausen  parameter  _A-  is  equal  to  -A.  * -1.36  x 10c 

2/5  , x 1/5  1/5 

This  Gan  also  be  written  as  _/\_  ck  rr  ( r ) R« 

U x 

J 

where  L * Lateral  integral  scale# 

11*  • X ^’t'2 

The  parameter  Tj”  1 r ) is  known  as  the  "turbulence  parameter". 

y ..t  ^ i/5 

Surprisingly  enough;  the  correlation  or  ( y — ) vs  . worked  very 

2/5 

well  for  the  existing  sphere  data  and  later  J\-  " was  used  to  correlate 

flat  plate  transition  data  by  Fage  and  Preston11  (Ra  transition  is  defined 


rrn#iiarNai)><rt>Mi 


- 30,  Re-  = 375 


10,  He-  = 67 


as  the  beginning  of  the  transition  region . } A = -12  corresponds  to  sep- 
aration in  the  Fchlhauscn  treatment;  the  subsonic  transition  data  give 
A - -2.  In  view  or  the  approximations  involved,  the  absolute  value  of 
A probably  doesn't  have  much  significance.  Figure  20  shove  A-A4"  vs . 
Re^tr  where  JX^'^oc  Hi.  ( “L.)-^  where  X^  is  the  longitudinal  microscale 

of  the  turbulence*  (For  the  present  problem  there  are  some  advantages 
in  using  X.  rather  than  h.)  In  addition  to  the  points  of  Hall  and 

12  ^ y 

Hislcp  (used  by  Fage  and  Preston),  two  values  from  Dryden’s  J work  are 
shown.  The  scale  of  the  turbulence  for  Dryden's  low  Reynolds  Number 
points  was  not  measured  and  was  not  easily  estimable  from  the  mesh  size 
of  the  screen  producing  the  turbulence.  A large  number  of  aluminum  flags 
were  tied  to  the  downstream  aide  of  the  screen  to  increase  the  turbulence 
level,  and  they  may  have  also  altered  the  turbulence  scale.  Because  of 
the  doubtful  value  of  X^  for  these  points,  the  curve  is  shown  dashed, 
Dryden’s  value  for  the  higher  Reynolds  Number  is  in  good  agreement  with 
the  values  of  Hall  and  Hislop. 

Dryden’s  value  was  the  lowest  transition  Reynolds  Number  point  this 
Author  was  able  tc  find  in  the  literature , It  is  naturally  interesting 
to  consider  whether  the  same  general  correlation  would  hold  for  still 
lower  transition  Reynolds  Numbers.  In  a restricted  sense,  the  laminar 
sub-layer  in  a fully  developed  turbulent  boundary  layer  can  be  used  to 
extend  the  correlation  curve  to  very  low  transition  Reynolds  Numbers, 


Laufer  4 has  shown  that  the  fluctuation  velocities,  as  well  as  the 

mean  velocities,  within  and  in  the  immediate  vicinity  of  the  laminar  sub - 

( yUr 

layer  when  plotted  as  vs  « — - are  independent  of  the  Reynolds  Number, 

0 V 

Since  the  mean  velocity  for  different  Reynolds  lumbers  follows  a single 
curve  of  jj-  vs  ~~  > if"  a value  for  is  decided  upon,  then  the  Reynolds 
Number  of  the  1 amine. r suo-Iciyei  bused  on  the  laminar  sub-layer  thickness  is 
invariant , i ,e . , 
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Tliis  laminar  sub-layer  can  be  thought  of  having  originated  at  soma  dis- 


. + . 6*1  '6 
tance  x upstream  wncrc  x - — -=r- 

vU-r  -J4 


S - y and  Ug  - U at-  — - 10  or  30  for  the  fiat  plate  la-dinar  boundary 


layer . Clearly  the  disturbance  in  the  turbulent  boundary  layer  us  l out- 
side of  the  laminar  sub-layer  is  so  intense  that  it  prevents  the  laminar 


sub-layer  from  growing,  that  is  in  terms  of  Reg . Tho  correlation  of 


( ^-)  vsc  Re  as  given  in  Figure  28  can  be  viewed  as  a purely 


empirical  result*  For  values  of  S-  ( ;■ — ) below  the  curve,  transition 


1/2 


does  not  occur.  The  value  of  this  3ame  parameter  can  be  computed,  fairly 
crudely,  for  the  laminar  sub -layer  uiing  tho  measured  proper+.i  ms  of  the 
boundary  layer  In  the  immediate  vicinity  of  the  laminar  sub-layer  as  the 
"free  stream" , The  parameter  can  be  only  loosely  determined  because  the 
boundary  layer  turbulence  is  anisotropic  and  there  are  many  fluctuation 
levels  and  many  microscales.  Taylor’s  value  for  the  spatial  pressure 
gradient  was  based  on  a Lagrangian  scale  Xrj  whose  relation  to  the  Eulerian 
microscale  was  determined  from  one  experiment  in  isotropic  turbulence. 
Intuitively,  the  most  important  determining  factor  in  a given  spatial 
gradient  would  be  the  appropriate  fluctuation  component  and  the  correspond- 
ing microscale,  such  as  u'  aid  Xx> 


yUr 


Using  u'  alone  and  Xx,  taking  the  edge  of  the  laminar  sub-layer  at 


- 30,  and  considering  the  "free  stream”  to  be  an.  average  of  conditions 


yU  yrr_*  ' 1/2 

from  • 30  to  — ^ « 60,  ( £ ) 


.31  Rextr  - li.l  x 10J. 


ytTr 


If  the  laminar  sub-layer  is  defined  as  extending  to 

similar  calculation  gives 

» 

( £—)  ■ .Up  where  Re. 


1C.  then  a. 


i /Vs  -*= 
X./C. 
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xtr 


>:  10' 
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The  values  of  for  the  laminar  sub-layer  are  used  to  extend  the 


correlation  curve  in  Figure  23  to  very  low  Reynolds  Numbers , 
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TURBULENCE  PARAMETERS  FOR  BOUNDARY  LAYER  Or!  CYLINDER 


In  order  to  compare  the  turbulence  parameter  =p  ( £ ) 


L/2 


ior  me  nw 


laminar  boundary  layer  along  the  cylinder  with  the  correlation  curve  in 
Figure  28,  the  turbulence  levels  and  microscales  for  the  possible  "free 
stream*  flows,  boundary  layer  or  wake,  are  required.  The  turbulence  level 
and  microscale  for  the  boundary  layer  free  stream  can  be  obtained  by  esti- 
mating the  change  in  turbulence  level  and  microscale  of  the  boundary  layer 
on  traversing  the  corner.  The  turbulence  level  and  microscale  for  the  wake 
have  to  be  based  on  low-speed  wake  measurements.  If  there  were  no  turbulence 
production,  then  the  boundary  layer  turbulence  at  the  beginning  of  the  cy- 
linder would  decay  rapidly  as  for  isotropic  turbulence.  Since  there  is 
turbulence  production,  the  turbulence  level  must  be  higher  than  that  obtained 
for  isotropic  decay.  On  the  other  hand,  the  turbulence  level  in  a wake  is 

u t 

limited,  and  asymptotically  -jj-  - 

8 * 


.2  where  U is  the  maximum  mean  3peed 
max 


decrement . 
decreases . 


Since  U^^  decreases  oc  x 
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, the  turbulence  level  steadily 
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will  be  C Clouted  fOi'-  the  following  cases. 

(The  Mach  Nurber  variation  along  the  cylinder  will  be  neglected.) 

(a)  The  mean  velocity  distribution  of  Figure  22  assumed  constant  and 
the  turbulence  assumed  to  decay  isotropically* 

(b)  The  mean  velocity  distribution  for  a wake  with  the  txrrbuler.ee 
level  and  microscale  based  on  low-speed  wake  data. 


TURBULENCE  LEVEL  AND  SCALE  IN  BOUNDARY  LAYER  AHEAD  OF  THE  CORNER 


Experimental  studies  of  low-speed  turbulent,  flows,  such  as  boundary 


-l  1 

X4 


layers, w channel  flow  and  pipe  flow,  **  have  established  the  following 
facts : 

(a)  Close  to  the  wall,  exit  to  the  region  where  the  velocity  distri- 
bution follows  the  logarithmic  law,  the-:  fluctuation  levels  for  all  cf  these 


flows  substantially  fall  on  a single  curve  when  the  data  are  plotted  in 

. u' 


terms  of  U-f,  th9  friction  velocity,  that  is  jj^.vs, 
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(b)  As  shown  by  Laufor  for  pipe  flaw,  well  away  from  the  wall  th® 
friction  velocity  is  the  proper  reference-  velocity  so  that,  a plot  of 

v s . ^ is  independent  of  Reynolds  dumber  r 
(c;  In  view  of  (a)  and  (b)  it  would  seen  reasonable  to  guess  that 
a curve  w^rs.  | for  the  boundary  layer  well  away  from  the  laminar  sub- 
layer would  be  odbsta  Hon+.  of  Reynolds  Number.  E.ebanoff  and 

i 

Diehl4  made  ilucoation  measurements  in  turbulent  boundary  layers  at  two 
different  air  speeds  and  at  several  different  boundary  layer  thicknesses. 

11  § 

While  all  of  the  data  plotted  in  terms  cf  do  not  fall  on  a single  curve, 
qualitative  results  are  consistent  with  the  hypothesis.  Furthermore,  near 
In  the  boundary  layer  well  away  from  the  wall  but  not  in  the  inter- 
mittently turbulent  region,  a coripar-5 son  of  KLebanoff'3^  boundary  layer  data 
with  Laufor’a^*  pipe  flow  data  in  the  same  region  shews  that  the  values  of 
^ are  approximately  the  sane.  For  these  computations,  the  fluouation 
levels  ^ as  measured  by  hlebanoff  are  used  for  the  portion  of  the  boundary 
layer  where  (a)  does  not  apply. 

(d)  The  experiments  of  hauler,  * Klebanoff  and  Diehl,4  and  Klsb&ncff , ' 
show  that  th?  longitudinal  integral  scale  of  turbulence  appears  to  be  con- 
trolled by  the  boundary  conditions  so  that  ior  a given  flow,  say  s boundary 
layer,  the  ratio  of  the  integral  scale  to  the  boundary  layer  thickness,  say 

Lx  ‘ Lx 

“g  is  not  a function  of  the  Reynolds  Number.  For  the  boundary  layer,  -y  - .4 


For  isotropic  turbulence  the  relation  between  the  microscale  of  tur- 

■'  6 

balance  and  the  integral  scale  is4 


7 


where  N« 


Since  the  turbulent  shear  flows  under  consideration  are  not  isotropic,  this 
relation  between  the  turbulence  scales  would  not  be  expected  to  hold . How- 
ever, a check  of  Klebanof f 1 s boundary  layer  data  shows  that  over  a fairly 
large  portion  of  tho  flow,  excluding  the  laminar  sub-layer,  the  isotropic 

relation  between  X a:,d  L holds  surprisingly  well » In  order  to  obtain 
- x 

suitable  values  for  the  turbulence  scales  in  these  high-speed  flows,  tho 

Xv  7 

approximation  that  >■  will  be  made. 

* 1 / % 


Hi 


At  the  cone  bass  at  M « 3.0k:,  with  H_  * 3.-3?,  Pn  m 130  cm.,  Ur**  2 .Ult.  x 10" 

_n  — ■ ^r* 

cm/sec.,  and  - 1 .1$  x 10  ' cm.  frigurs  29  shows  the  turbulence  levels  tj  , 
across  the  boundary  la"?*'  "V^rn  n is  the  mean  velocity  at  ary  point  across  the 
boundary  layer  and  u*  is  the  square  root  of  the  average  of  the  squares  of  the 

i i i 

turbulent  velocity  components  u , v , and  v . 

THE  EFFECT  OF  THE  EXPANSION  AT  THE  SHOULDER 
ON  THE  BOUNDARY  LAYER  TURBULENCE 

Consider  a flew  in  which  there  is  a vorticity  distribution  super- 
imposed on  a uniform  mean  motion.  When  there  is  an  irrotational  change  of 
the  mean  motion,  tne  rotational  and  irretc.tional  flows  will  interact,  so 
that,  in  general,  the  energy  of  t!re  rotational  field  is  changed.  In  con- 

17 

flection  with  the  design  of  wind  tunnel  contractions,  Prandtl  • showed  that 
if  the  gain  in  energy  for  all  of  the  streamlines  in  a contraction  is  the 
same,  a condition  that  would  realized  if  the  pressure  drop  for  each  stream- 
line were  the  same,  then  a longitudinal  steady  perturbation  in  velocity 
would  be  reduced.  He  also  applied  the  circulation  theorem  to  a vortex  with 

its  axis  in  the  stream  direction  and  showed  that  its  energy  would  be  in- 

1 8 

creased  in  traversing  the  contraction.  G.  I.  Taylcr,  and  more  recently, 

19  20 

Ribner  and  Tucker,  and  Tucker  have  investigated  the  interaction  of  an 
Irrotational  and  a rotational  flow  field  using  the  Helmholtz  Equations, 

The  Helmholtz  Equations  give  the  vorticity  of  a fluid  element  in  terms  of 
its  vorticity  at  some  initial  time  and  the  motion  of  the  fluid  element  with 
respect  to  its  initial  coordinates,  and  so  it  i3  a Lagrangian  representation. 
The  modification  of  the  vorticity  distribution  by  viscosity  is  neglected  so 
that  viscous  diffusion  and  a<?ous  dissipation  are  not  included.  The  equations 
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where  £q,  are  the  initial  components  of  the 

vorti-itv  vcc‘t>rYr'' 

a,  b,  c,  are  the  initial  cartesian  coordinates 
of  the  fluid  element 
pA  *•  initial  fluid  element  density 

v 

and  they  are  valid  for  a compressible  fluid.  In  a general  case,  to  obtain 
the  vorticity  components  at  time  t,  all  of  the  nine  Lagrangian  coordinate 
derivates  have  to  be  known . The  above  authors  simplify  the  problem  by 
assuming  that  the  changes  in  the  potential  motion  are  large  and  rapid,  and 
the  rotational  motion  is  relatively  small  30  that  the  motion  of  a fluid 
element  is  primarily  specified  by  the  velocity  changes  of  the  potential 
motion*  Under  tbit  restriction,  six  of  the  Lagrangian  derivatives  become 
aero  and  only  , and  remain,  sc  that  the  equations  describe  the 

alteration  of  a vorticity  distribution  by  a geometrical  deformation  of  the 
flow  field*  Ary  desired  vorticity  distribution  can  be  represented  by  the 
superposition  of  elementary  vorticity  distributions. 


These  methods  have  been  developed  to  deal  with  the  problem  of  the 
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changes  of  turbulence  in  a wind  tunnel  contraction  where  unfortunately 

in  general,  all  of  the  Lagrangian  coordinate  derivatives  are  important 

* 

and  where  the  dissipation  cannot  be  neglected.  As  previously  noted,  in 
the  present  experiments,  for  a good  portion  of  the  boundary  layer,  the 
acceleration  at  the  shoulder  of  the  cone  cylinder  model  occurs  in  a length 
of  the  order  of  6 so  that  compared  to  the  usual  wind  tunnel  contraction 
the  changes  in  the  mean  motion  of  the  turbulent  bewndary  l^yer  fluid  occur 
very  rapidly.  The  following  estimates  for  the  importance  of  the  non-linear 
terms  and  the  viscous  dissipation  in  the  expansion  at  the  cone  cylindor 
• shoulder  indicate  that,  in  this  case,  the  linear  theory  should  be  reasonably 
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For  the  turbulent  boundary  layer  at  the  base  of  the  cone  L^  - .US  and 


«•  = .0?  (Average  from  % = .1  ^ *•  „U) . 
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An  average  value  of  ^ for  the  same  region  is  jt- 


See  Appendix 
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transverse  to  the  direction  of  the  mean  motion,  then  the  motion  of  a fluid 
element  depends  on  its  initial  y coordinate  b,  as  follows: 

x - U x At  + j/a  + ( TT"'1  x b 

mean  t U mean 


TT  mean  velocity  in 
mean 

acceleration 


, Sx  « dx  6 f vl'  \ 

81x1  35  ’ * ’ 35  * IT  ( lU 
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3x  /ox 

3b  /5a 


, a 2 aver  ^ , r 3x  /ox  ^ 9 

where  f 1 ~ l.->  or  15b  / cia  ~ 

1 aver  ' 

so  that  even  for  this  case,  the  neglected  derivatives  are  not  negligible 
although  they  arc  smaller  than  the  derivatives  that  are  included,  and  the 
assumptions  in  the  calculations  are  at  least  roughly  satisfied.  It  is  also 
necessary  that  the  turbulent  dissipation  3.n  the  flow  about  the  shoulder  be 
unimportant  as  compared  to  the  change  in  turbulent  energy  produced  by  the 
fluid  element  defoliation. 


The  decay  of  the  turbulence,  if  isotropic,  would  be  proportional  to 

,2  ~«3 


d(u  ) 
dt 
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The  time  for  passage  past  the  shoulder,  assuming  the  acceleration  length 
is  6 is  At"  — 


Then  the  loss  of  turbulent  kinetic  energy  is  (AuV  = - ^ 


mean  6 


The  fraction  of  the  turbulent  energy  lost  is 


&u'  )2 
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„ mean  7 

Kcw  —a*  « .25  so  that  — ■>. — due  to  fluid  element  deformation  is  ,75  as 

~»2  u«2 

1 1 

compared  with  .06  for  viscous  dissipation. 

Uk 


I 


Ti  has  been  shown  that  the  restrictions  of  the  linearized  theory 
sre  fairly  well  satisfied  for  ih*  rapid  acceleration  at  the  shoulder 
of  the  cone  cylinder.  Since  the  available  measurements  on  the  overall 
changes  in  turbulence  intensity  in  wind  tunnel  contract!  <-r.m  are  in 
general  agreement  with  the  calculations,  it  is  believed  Justifiable  to 
use  the  theory  for  obtaining  the  overall  turbulence  intensity  changes 
at  the  shoulder.  However,  certain  inaccuracies  are  manifest*  The 
turbulent  field  ahead  of  the  corner  is  neither  isotropic  as  in  reference 
19  nor  axisymmetric  as  in  reference  20.  Furthermore,  the  expansion 
right  at  the  model  shoulder  is  two-dimensional  rather  then  axisymmetric 
which  is  the  case  computed  in  the  above  references. 


The  conclusion  as  to  the  general  nature  of  the  phenomenon,  on  the 
cylinder  is  not  dependent  on  on  accurate  knowledge  of  the  turbulence 
level.  Since  many  other  approximations  are  made,  it  did  not  seem 
worthwhile  to  attempt  to  evaluate  the  possible  error  due  to  the  above 
differences.  It  is  noted  by  Ribner  and  Tucker  that  their  result  is 
very  el-036  to  the  result  obtained  by  Prandtl  for  a much  simplified 
vorticity  distribution.  In  any  case,  the  author  could  not  resist  the 
pleasure  of  merely  plucking  an  attenuation  factor  from  a handy  chart. 
Figure  29  shows  the  fluctuation  levels  after  the  expansion  at  the 
shoulder  as  obtained  using  reference  19. 


APPLICATION  TC  TURBULENCE 

figure  30  shows  an  expansion  of  a turbulent  flow  where  the  change  in 
mean  velocity  is  large  and  rapid  enough  to  permit  the  use  of  Helmholtz's 
Equations  with  G.  I.  Taylor's  restrictions.  The  r.r,T cnlaticn  can  be  described 
in  the  following  terms . 


At  any  time  t,  an  instantaneous  measurement  of  the  velocity  field,  if 
it  could  be  made _ would  establish  the  vortioity  of  the  fluid  elements  at 
station  A.  The  Helmholtz  Equations  can  Chen  oe  used  to  determine  the 
change  in  vortic ity  of  the  fluid  elements  at  station  B.  If  a sufficient 
number  of  similar  measurements  are  made  randomly,  the  average  changes  in 
vorticiiy  of  the  fluid  elements  can  be  obtained,  and  so  the  average  change 
in  energy  in  the  relational  motion. 


The  turbulent  energy  equation  represents  the  change  in  turbulent 
kinetic  energy  in  a fluid  element  moving  with  the  mean  motion.  Ary  par- 
ticular fluid  element  may  lose  or  gain  turbulent  energy  in  various  ways, 
such  as  turbulent  energy  convection,  but  the  fluctuation  field  as  a whole 
can  only  have  its  total  turbulent  kinetic  energy  changed  by  losing  kinetic 


energy  to  heal  through  viscous  dissipation  or  by  having  its  kinetic  energy 
changed  by  interaction  with  the  mean  motion. 


The  interaction  term,  which  is  usually  called  the  turbulent-  production 

au . 

term  is  -pu^  1 


*» 


(1,  j are  repeating  indices) and  may  be  either 


positive  or  negative,  ’'/hen  the  mean  motion  is  uniform,  then  the  production 
term  is  zero.  If  the  viscous  dissipation  is  negligible  in  the  section  A-;3, 
then  the  integration  of  the  production  term  over  the  volume  A-B  must  result 
in  a change  in  turbulent  kinetic  energy  equal  to  the  change  in  turbulent 
energy  calculated  using  the  Helmholtz  Equations. 


Ths  calculations  of  Ribner  and.  Tucker-1"7  show  that  mean  flow  changes  may 
result  in  either  a decrease  or  an  Increase  of  turbulent  energy*  For  acceler- 
ations at  supersonic  speeds,  the  fluid  elements  are  stretched  in  both  the 
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loogjtudl nal and  the  lateral  directions.  As  a result,  both  the  longitudinal 
and  lateral  fluctuation?  arc  reduced  by  the  defamation  and  the  turbulent 
energy  is  thereby  reduced.  In  a senna,  the  turbulent  energy  flows  fres  the 
tuj^nlAnt  fli»rttTtation  f iftld  to  th2  l»"»"  floyt  This  conversion  Of  rfUadOC 
-""-r”.  *— »<>  ^nalogcus  to  ths  conversion  of  rsnacw  ^henu.'i 

. i 

energy:  in  a gas  into  directed  motion  in  a supersonic  nozzle.  If  the  flow 
process  is  reversed,  the  turbulent  energy  will  be  increased* 

' STAGNATION  TEMPERATURE  UNEVENNESS  IN  THE  TURBULENT  BOUNDARY  LAYER 

Kcvassnay5^  hot-wire  measurements  in  a supersonic  boundary  layer  at 
M ■ 1*75  show  that  adjacent  fluid  nasces  vary  in  stagnation  temperature! 
if  a measurement  is  made  at  a fixed  point  in  the  boundary  layer  there  is 
a stagnation  temperature  fluctuation.  The  magnitude  of  the  fluctuation 
varies  through  the  boundary  layer  end  reaches  a peak  value  of  about  % at 
the  gtagr.9-M.cn  temperature , As  Kovasziviy  discusses  in  his  paper,  the  stat- 
nation  temperature  fluctuation  involves  a density  fluctuation  but  does  not 
involvs  a pressure  fluctuation. 

However,  when  a field  of  uneven  stagnation  temperature  is  accflerat-ed 

to  higher  speed,  a vorticlty  field  consisting  of  longitudinal  velocity  v ari- 

?3 

atioss  should  b9  created.  C cursin'"  has  derived  an  expression  for  the 
magnitude  of  the  resul ting  velocity  perturbation  in  terms  of  the  initial 
stagnation  temperature  perturbation.  The  relation  is 
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where  • stagnation  temperature  perturbaticn  * 


If  Kovasznay's  measurements  are  applied  to  the  turbulent  boundary 
layer  «t  “ 1.8?,  we  have,  as  maiiciimaa  velum*  of  tlie  root  mean  square 
fluctuations 
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Assuming  that  this  specially  created  vorticity  field  would  become  a part 
of  the  jecaral  vorticity  field  after  the  model  shoulder,  the  contribution 
to  the  turbulent  fluctuation  level  would  be  .00b . This  addition  tc  the 
turbulence  after  the  shoulder  will  be  neglected. 


EFFECT  OF  SHOUIDER  EXPANSION  OH  TURBULENCE  HICRCSCALE 

For  isotropic  turbulence,  the  integral  scales  are  often  considered 

20 

to  define  a sort  of  average  eddy  size.  Tucker  has  calculated  the  eifect 
of  stream  velocity  changes  cn  the  correlation  i icxcn m h um  wmen  the 

integral  scale  is  determined.  He  notes  that  in  certain  cases  the  longi- 
tudinal scale  becomes  negative  which  .nakes  fairly  dubious  the  identification 
of  the  integral  scale  as  an  average  eddy  sise. 


The  turbulence  microscale  can  be  defined  as  X 


2 _ 2u» 
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so  that  the 


change  in  X is  essentially  determined  by  the  stretching  in  -one  x direction. 


that  is 


X x 1.U2 


1 refers  to  position,  just  ahead 


2 refers  to  position  ju3t  aftor 


corner 


But  in  the  expansion  at  the  shoulder,  the  expansion  of  tho  fluid  elements 
in  the  y direction  is  not  equal  to  that  in  une  x direction,  so  that  with 
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At  the  beginning  of  the  cylinder  tee  "free  stream”  turbulence  in  these 
experiments  is  certainly  anisotropic)  even  if  it  were  isotropic  ahead  of 
the  shoulder  it  would  bo  anisotropic  after  the  expansion.  The  flow  near 
the  center  of  a two-dimensional  wake  is  roughly  isotropic,  so-  thot  a*  the 
"free  stream"  flew  takes  on  the  aspect  of  a wake,  field 
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3hould  move  towards  isotropy.  Since,  in  ary  case,  the  turbulence  in  the 
’’free  stream"  along  the  cylinder  will  be  compared  with  low-speed  data  where 
the  free  stream  turbulence  vas  isotropic,  the  anisotropic  turbulent  field 
at  the  beginning  of  the  cylinder  will  be  described  in  terms  of  an  equivalent 
isotropic  flow.  Hence  the  geunetrlo-1  charge*  in  the  expansion  will  be 
averaged,  giving  X„  ■ 1.5U  X,  . 
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BOUNDARY  LAYER  "FREE  STREAM” 
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Following  are  average  values  for  the  region  in  the  boundary  layer  ahead 
of  the  corner,  from  ^ - .2  to  ? - Jj. 

X'  “ 3-3  x ICp  cm/sec  X - 1.15  x 10~2  cm. 

X 

x 

Altar  the  corner, 

ul  - 2.1  x 10^  cm/sec  X - 1.76  x 10~2  cm. 

-2  *2 

To  estimate  the  turbulence  decay  in  the  "free  stream"  along  the  cylinder, 
it  is  necessary  to  establish  the  integral  seals  or  an  effective  mesh  ?i«&  for 
the  equivalent  isotropic  flow.  This  can  be  dons  as  follows. 


For  isotropic  turbulence,  the  decay  of  the  turbulent  energy  is  given  by 
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which  is  based  on  the  expression  for  the  viscous  dis- 
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sip  at  ion  W,  where  V - 

X 41 


Except  near  the  laminar  sub-layer,  the  dissipation  has  been  found  to  be 
reasonably  well  represented  by  the  above  formula.  Now,  from  measurements 
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behind  grids  of  mesh  size  Meff’  i+'  has  been  found  that1 
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which  can  be  written  as 
* Meff  US 
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An  effective  mesh  size  for  one  equivalent  isotropic  flow  can  be 
determined  by  equating  the  above  expressions  for  the  nurbuience  decay . 


Then  after  the  corner  * 5 S * 1 O cm.  and  it  follows  that 
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The  results  of  the  calculations  for  the  turbulence  level  and  the 
microscale  are  shewn  in  Figure  31.  The  resulting  turbulence  parameter 

u/v  V2 

jy  f 57”“)  sll0wn  i*1  Figure  32.  At  any  position  along  the  cylinder, 

' x' 

x,  as  used  here,  is  the  affective  length  of  the  new  laminar  boundary 
layer  at  that  position,  so  that  is  approximately  equal  to  the  distant- 
from  the  beginning  of  the  cylinder. 


WAKE  FREE  STREAM 


The  only  measurements  or  trie  i.iu-b./Leriee  in  two-dimensional  wakes  that 

g 

the  author  wa3  able  to  find  are  those  of  Towsend’"’  where  the  distributions 

of  u*  and  the  microscale  \y  were  measured.  Previously  it  was  noted  that 

for  such  flews  as  channel  flows  and  boundary  layers,  the  integral  scale  L 

L.  _ x 

deta-r«-<-5 ned  bv  the  bound arv  conditions,  so  that  — was  a constant  in- 

6,, 

dependent  of  Re.  It  will  be  assumed  that  for  a wake,"  also,  the  integral 
scale  is  a constant  fraction  of  the  wake  width,  so  that  we  assume  — constant. 

yo 

From  Townsend,  for  Uc  *■  1.12  x 10J  cm/sec. 

-1 

t~“  - 1.30  ax*l  -2  = 3.18  x 10  x 
x y0 

(y^  = half  wake  for  mean  velocity  distribution) 
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So  — ■ = .2l4i 

y0 

Since  the  equivalent  wake  at  the  beginning  of  the  cylinder  has 
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X ■ 3.f>0  x 7l0-‘  cm. 
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•^Suggested  by  Professor  Stanley  Gorrsin  of  Johns  Hopkins  University 
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which  compares  with  ■ 1 ,?6  x 10~*  cm.  for  tha  boundary  layer  after 
the  corner. 


The  decay  as  the  wake  width  grows  can  be  obtained  vising  ths  fact 

tha*  ^222  « i .*12. 
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where  x ■ Xe.  ■*•  x* 

w «» 

x_  Effective  starting  length  for  the  wake 
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u*  x 2y  a:T  distance  -rorn  corner 
Since  Rst  - — — ■ constant  for  ary  given  wake, 

L u*Lv 

and  — * constant,  then  — - — ■ constant. 
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or  Xx  - .15  yQ 
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Also,  we  havs  y^  • c^x  , o.  « yonetant. 


Tha  turbulence  levels  and  micros oals  values  are  shown  in  Figuz-e  31* 

The  resulting  "turbulence  parameter"  is  shown  in  Figure  32 • Allowance  has 
been  made  for  the  loss  of  wake  fluid  to  a new  laminar  boundary  layer. 

Ths  start  of  the  abrupt  aurfaoe  temperature  rise  on  the  cylinder  is 
not  a well  defined  point  sines  the  measuring  points  were  not  closely  spaced* 
The  position  marked  on  Figure  32  as  representing  transition  on  the  oyiinder 
is  2*5  cm.  from  the  corner* 


In  a low-speed  wind  tunnel  test  the  frse  stream  turbulence  is  of  rela- 
tively large  scale  so  that  ever  the  length  of  a flat  plate  model  in  the  teat 
section  the  changes  of  the  turbulence  level  and  turbulence  scale  ftps  small* 
The  tusbulencs  parameter  I"  ur  uhjG  lsv---speed  boundary  layer  incraasss  steadily 
proportional  to  and  would  necessarily  intersoot  the  oorrslation  curve* 
In  tha  present  case,  ths  turbulence  soale  in  the  "free  stream"  along  the 
cylinder  is  relatively  small  so  that  significant  decay  ef  the  turbulence 
occurs  il one  tha  run  of  tha  new  laminar  boundary  layer.  Apparently,  the 
calculated  deosy  is  large  enough  so  that  the  turbulence  par erne ter  for  the 
boundary  layer  free  stream,  decaying  isotropically,  has  a maximum  value 
along  tho  cylinder.  Figure  32  shows  the  turbulence  parameter  curve  just 
touching  the  correlation  curve. 


Actually,  as  has  been  noted  previously,  the  turbulence  level  of  the 
"free  stream"  should  decrease,  mere  slowly  than  for  isotropic  decay  because 
of  the  turbulence  production  taking  place  in  the  "free  stream"  mixing 
processes-  At  the  beginning  of  lh«  cylinder  the  turctuen««  level  In  the 
boundary  layer  "free  stream"  is  higher  than  that  for  the  equivalent  wake 
"free  stream"  at  the  same  position  as  shown  in  Figure  31.  bVsniually,  the 
turbulence  parameter  for  the  wake,  will  also  reach  a maximum  and  then  de- 
crease, hut  ovsr  the  region  of  interest  it  increases  monotonically.  it 
appears  likely  that  the  turbulence  parameter  for  the  boundary  layer  on 
the  cylinder  is  either  on  or  above  the  turbulence  parameter  for  the  wake. 

But,  in  ary  case,  it  is  clear  from  the  curves  in  Figure  32  that  the  general 
nature  of  the  turbulence  parameter  along  the  cylinder  is  the  3ame  for  both 
of  the  "free  stream"  assumptions.  The  confluence  of  all  of  the  curve;?  in 
the  vicinity  of  the  experimental  value  of  Revtr  would  appear  tc  support  the 
hypothesis  concerning  a laminar  turbulent  transition  on  the  cylinder. 

Figures  16  and  17  show  that  there  arc  differences  in  the  temperature 

u*  1 

data  at  different  tunnel  pressure  levels.  If  the  value  of  x — y were 

A 

x 

the  same  at  each  pressure  Isvsl  | v hen  x for  uhc  beginning  of  transition  on 
the  cylinder  would  be  the  same  at  all  of  the  pressure  levels,  and  Re^.^ 
would  increase  proportional  to  the  pressure  level.  At  M ■ 3*02,  this  is 
roughly  the  experimental  result.  At  M ■ 3«5>£  only  the  100  cm.  data  has 
transition  oeginning  at  a significantly  different  position  along  the  cylinder 
as  compared  with  tho  data  at  the  other  pressure  lsvals . The  sizes  of  the 
wire  trips  used  on  the  cone  at  the  different  pressure  levels  were  not  care- 
fully controlled  to  make  the  boundary  layer  conditions  at  the  trips  "similar1" 

u1  1 

as  the  pressure  level  was  varied.  Since  the  narameter  x — will  de- 

x 

pend  on  both  the  turbulent  boundary  layer  thickness  at  the  cone  base,  and  the 

tunnel  pressure  level.  x — i,— would  be  expected  to  vary  in  a random 

' u X 

x 

way  with  tunnel  pressure  level. 


EFFECT  OF  CKlITtESSIBILITY 


Ociftpressibiliiy  effects  rave  ueen  included  for  that  portion  of  the 
anajysia  involving  the  mean  flew  properties  ahs<wi  of  the  comer  and  on 
the  cylinder;  The  introduction  of  turbulence  levels  and  turbulence  0 coles 
was  based  on  low-speed  data.  For  this  reason  alone,  the  close  agreement* 
between  the  predicted  value  of  Re  . and  the  experimental  value  of 

XT»r 

could  be  considered  fortuitous. 


T'ns  supersonic  boundary  layer  measurements  of  Kovasziwyf  can  be  com- 
pared with  the  low-speed  fluctuation  level  measurements.  Kc vasznsy  made 
his  measurements  at  the  following  conditions 
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So  the  value  of  n»  expected  on  the  basis  of  the  low-speed  measurement 
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temperature  should  be  below  room  temperature,  so  that  the  wall  temperature 
will  be  raised  somewhat  by  heat  transfer  from  the  room*  A higher  wall 
temperature  would  increase  U-£  and  thereby  u’ ) • 

The  influence  of  Mach  Number  on  the  "turbulence  parameter"  can  be 
obtained  by  making  use  of  an  existing  solution  for  'nigh-speed  boundary 
layers.  Doridnitsen.^  under  the  restrictions  that  the  Prandtl  Number 
eoualed  one  and  that  the  heat  transfer  is  zero,  developed  a method  for 
handling  a compressible  boundary  layer-  analogous  to  the  Karman-rohlhauson 
method  for  a low-speed  boundary  layer. 


ihe  boundary  layer  profile  is  expressed  as  a _e 

W * AT  + Bt2  + Ct3  + DT^* 
where  A ■*  ( 2 * ) , B - - 
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V ~ — where  i is  a transformed  coordinate  of  the  distance  from  the  wall  y 


as  follows 


dt  - r - dy 


- w(  -c) 


L - characteristic  length  cf  problem 
6 = edge  of  boundary  layer  in  coordinate  t , 
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= 0 gives  _/L  - - 12  for  separation  as  for  the  low-speed 


solution = 


The  parameter  S\~  is  as  fellows t 


_A_  - ZSL  2 -j~  j Vg  - — c - velocity  of  efflux  into  vacuum 
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s is  a transformed  x coordinate,  where 


<«  local  free  stream  pressure 
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For  a compressible  fluid 
usdUs  ♦ 0 


' s ~s 

Sin<:e  p0  ns  po 

^ . i.  ^ dt  , 6-i  J ^ * 


or  with 


dt  ■ 6 d T 


.«/§ dr > F-^v] 


" “ Y 

r f->  *.  Y“I  «2N,  Y’  C f-\  -1>  nr  Ur 

" VJ  T‘V  o7  V"  » V 


SuDStitstisg  for  6 and  vg  nave 

ds 


_/U  vfe  h(%>  A'7  ($rl 


| 


3 

\ 

I 


s 

<■« 


r 


. nl  2,1-2 

-L<y  a - Us*1)  dTi 


"c"  denotes  compressible 


xms  can  oe_comparea  wion 


_5<:  / axt  \ 

— “ j for  an  incompressible  fluid 


inc  n, 


* inc *'  udiotoo  iuc w» awCibl 


U.-!) 


so  _/L„  “►A 

r$ 


Now  if  a turbulence  field  is  carried  along  by  a supersonic  flow.  them 
should  be  instantaneous  pressure  gradients  analogous  to  those  for  a low-speed 
turbulence  field. 


We  then  have,  simply  assuming  that  the  low -speed  relation  holds 


lio„u1'1 

•j 


A -6^ 

or  A„  - 7-7?“  x 

c uobs 


For  a high-speed  boundary  layer 
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If  it  is  assumed  that  the  same  value  of  j\.  , compressible  ox-  incom-  , 

u*  / x v' ^ 

nressible.  is  required  for  transition,  then  a smaller  value  of  JJ~~  ( J~) 


o \ x/ 


is  required  for  the  compressible  fluid  since  .fi(l-L)  is  in  general  gj.AAt.pw 

then  uni  tv.-.  Figure  33-  To  evaluate  G,  , was  varied  from  G to  **12.) 

So  Re  , at  high  speeds  will  be  less  than  at  low  speeds  for  the  same  value 
xtr 
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for  transition.  Rut-  Kov-sz-nay’s  results  show 
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of  ,7»  so  that  the  two  different  effects  of  compressibility  would  approxi- 
mateiy  cancel  in  this  experiment. 


If  the  concept  of  transition  being  caused  by  the  instantaneous  pressure 
gradients  due  to  ulv©  Mui’l  julent  field  has  aiy  validity,  then  the  superposition 
of  a steady  pressure  gradient  of  sufficient  magnitude  shauld  prevent  the 
occurrence  of  instantaneous  adverse  gradients . In  Figure  32  the  calculated 

‘ bb-l...iCo  par Aaw  for  the  boundary  layer-  flow  from  the  beginning  of  the 

cylinder  follow  fairly  close  and  almost  parallel  to  the  correlation  curve. 

The  decay  of  the  turbulence  in  the  "free  stream"  is  very  significant  and 
the  turbulence  parameters  for  the  new  laminar  boundary  layer  tend  to  level 
out  fairly  quickly.  The  present  model  actually  has  a 3inall  adverse  pressure 
gradient  which  has  been  neglected  in  these  calculations.  It  seems  possible 
that  a model  designed  with  a suitable  favorable  pressure  gradient  after  the 
shoulder  might  have  much  longer  runs  of  new  laminar  flci-j  since  the  favorable 
p2?Cw>s‘ux*s  gr^dioriv  rnight  X 0*702?  *bhs  \irh ilL 0 n c 0 p3.i,?jn6t  3 t cih*vq  for  "t^-3  boundfurv* 
layer  enough  to  avoid  transition. 


Nothing  in  the  present  analysis  indicates  that  the  pressure  drop  need 
be  a particularly  sudden  one*  In  a gradual  turn,  the  laminar  boundary  layer 
would  grow  and  the  turbulence  level  would  be  steadily  reduced.  It  5s  neces- 
sary that  the  acceleration  keep  the  turbulence  parameter  at  ary  point  at  a 
respectable  level.  Tests  on  a modsl  with  a .3  cm.  R at  the  shoulder  gave 
the  same  results  as  for  a 3harp  cornered  model.  Spreading  the  acceleration 
over  a much  greater  length,  perhaps  up  to  one  or  two  centimeters,  may  still 
produce  a new  laminar  boundary  layer. 


s 


i. 

£ 


CTf-fi'**  uv 


In  this  work  it  has  been  necessary  to  make  a distinction  between  the 
boundary  iaytir  a l a particular  position  on  a model  and  the  complete  shear 
layer  at  the  same  position*  The  complete  shear  layer  represents  the  results 
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the  model  surface  from  the  leading  edge  or  tip.  The  boundary  layer > as  it 
is  used  here,  is  defined  in  a local  sense,  and  may  be  all  or  only  the  inner 
portion  of  the  complete  shear  layer.  In  the  latter  case,  the  outer  portion 
of  the  shear  layer  is  the  outer  flew  or  environment  in  which  the  boumaxy 
layer  grows . As  such,  it  may  have  substantial  influence  on  the  friction 
a.id  neat  transfer  as  well  as  the  state  of  the  boundary  layer. 


A laminar  boundary  layer  is  considered  present  when  the  knowledge  of 
the  surface  friction  at  a particular  position  i?  being  transmitted  to  the 
outer  flow,  which  may  be  the  outer  part  of  the  complete  shear  layer,  by 
molecular  viscous  forces.  For  instance,  in  Figure  3k,  the  velocity  change 
of  the  wall  is  propagated  out  into  the  flow  by  molecular  viscous  forces  and 
so  the  boundary  layer  at  that  position  is  laminar. 


A turbulent  boundary  layer  exists  when  the  propagation  of  the  surface 
friction  at  a particular  position  on  the  surface  to  the  outer  flow  is  con- 
trolled by  turbulent  shearing  stresses.  Of  course,  even  when  the  boundary 
layer  is  turbulent,  there  is  a laminar  sub-layer  adjacent  to  the  wall,  but 
some  information  about  the  3 hear  at  the  wall  is  very  rapidly  communicated  to 
the  whole  turbulent  boundary  layer.  The  growth  of  the  turbulent  boundary 
layer  at  any  position  can  be  described  in  terms  of  th6  wall  friction  at  that 
position,  so  that  the  turbulent  boundary  layer  behaves  an  an  entity. 


Hence,  the  label  laminar  or  turbulent  is  selected  on  the  basis  of  how 
the  state  of  shear  on  the  surface  is  being  propagated  or  transmitted,  away 
from  the  surface  to  the  outer  flow. 
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The  growth  of  a boundary  layer  cn  a flat  plate  in  a wind  tunnel  air- 
stream  may  alter  the  properties  of  the  airstream  in  several  ways . Besides 
potential  changes  of  the  free  stream  which  are  dependent  ™_  the  bcunciarT 
layer  displacement  thickness  * it  is  also  possible  for  lias  boundary  layer 
to  have  a more  subtle  influence  on  the  free  stream  flow*  As  the  boundary 
layer  grows,  turbulent  air  is  being  removed  from  the  main  stream  and  either 
has  its  turbulent  energy  dissipated  when  the  plate  boundary  l^yer  is  laminar, 
or,  in  the  usual  case,  has  it  increased  if  the  plate  boundary  layer  is  tur- 
bulent.  In  either  case,  the  removal  of  a layer  of  turbulent  air  from  the 
main  stream  may  affect  the  main  stream  in  the  following  way* 

The  wind  tunnel  air  is  turbulent,  and  if  it  is  also  non-uniform,  then 
turbulent  mixing  will  cause  a progressive  alteration  of  the  mean  velocity- 
distribution  and  the  turbulence  of  the  air  as  it  passes  tlirough  the  test 
section.  The  removal  of  a layer  of  test  section  air  must  necessarily  change 
the  course  of  the  turbulent  mixing  processes . In  the  usual  wind  tunnel  case, 
these  affects  are  negligible  since  a serious  attempt  is  usually  made  to 
achieve  uniform  flow  in  the  test  section.  In  the  flow  about  the  cone-cylinder 
model  used  in  these  experiments,  where  the  outer  portion  of  the  complete  3hear 
layer  at  a particular  position  constitutes  the  "free  stream"  for  the  boundary 
layer,  these  effects  have  had  to  be  considered. 

Once  the  boundary  layer  at  a particular  positron  has  buen  identified 
as  being  laminar  or  turbulent  in  the  local  sense  defined  above,  then  the 
determination  of  the  local  friction,  and  heat  transfer,  requires  s knowledge 
of  the  history  and  development  of  the  entire  shear  layer  near  the  model 
surface • 

In  this  report  evidence  has  been  presented  to  =;how  that  along  a model 
surface  the  boundary  layer  may  be  successively  laminar,  turbulent,  laminar, 
and  turbulent  again.  It  was  found  that  with  a laminar  boundary  layer  at 
the  base  of  the  cone,  the  boundary  layer  on  the  beginning  of  the  cylinder 
was  also  laminar,  but  that  the  strong  pressure  drop  at  the  shoulder  had  a 
lasting  effect  on  the  subsequent  boundary  layer  development.  For  a rela- 
tively large  distance  downstream  of  the  pressure  drop,  it  was  not  sufficient 


to  describe  the  boundary  liver  properties  s n tews  of  the  local  '-'ash  Number 
just  outside  of  the  boundary  layer,  but  rather,  tha  history  of  the  flow  had 
to  be  considered. 

When  the  boundary  layer  was  turbulent  at  the  cone  base,  the  result  cf 
the  strong  pressure  drop  was  to  allow  the  growth  of  a laminar  boundary  layer- 
«•  t/  whu  bC/j  inning  of  the  cylinder.  This  new  laminar  boundary  layer 
subsequently  underwent  transition  to  a turbulent  boundary  1 ay or  a short,  dis- 
tance back  on  the  cylinder.  The  recovery  factors  along  the  cylinder  wero 
found  to  depend  on  the  local  laminar  or  turbulent  state  of  the  boundary  layer 
and  also  upon  the  history  and  development  of  the  shear  layer  near  the  model 
surf  ace  , These  two  cases  are  portraved  in  Fi  gure  3v  -• 

The  second  tradition  from  laminar  to  turbul  ant  flow  which  occurs  on 
the  «yl  liiier  of  the  cone -cylinder  model  i3  believed  to  be  controlled  by  the 
high  turbulence  level  of  the  " free  stream*'  for  the  new  laminar  boundary 
layer  starting  at  the  beginning  of  the  cylinder.  The  experimental  data  are 
compared  with  existing  lev -speed  correlations  for  laminar-turbulent  tran- 
sition in  highly  turbulent  free  streams.  Further,  a compressible  analogue 
for  Taylor's  turbulence  parameter  is  derived  and  it  is  found  that  at  least 
for  the  present  experiments  the  error  in  using  low-speed  turbulent  boundary 
layer  data  and  the  low-speed  correlation  data  should  not  be  large. 
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While  these  experiments  are  at  supersonic  speeds,  the  same  phenomena 
would  be  expected  to  occur  at  low  speeds  given  suitable  conditions.  For 
instance,  wind  tunnel  contractions  have  a strong  pressure  drop  and  as  a 
result,  although  the  boundary  layer*  in  the  settling  section  is  turbulent, 
the  boundary  layer  on  the  nozzle  wall  may  be  laminar. 
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corner  expansion  on  the  temperature  recovery  factors  were  made  in  associ- 
ation with  Mr.  B.  des  Clers . Recovery  temperature  variations  qualitatively 
similar  to  those  shown  herein  were  obtained,  but  the  situation  was  confused 
by  th9  large  heat  transfer  in  the  model  resulting  from  its  construction, 
and  from  our  inability,  at  that  time,  to  satisfactorily  reproduce  the  re- 
cmroTy  factor  levels  on  successive  tests.  The  lucite  model  was  designed  to 
greatly  lessen  the  model  heat  transfer. 


The  author  wishes  to  thank  Miss  Carol  Sande  for  her  substantial  help  in 
transcribing  the  manuscript  and  in  proofreading  the  report . 
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APPENDIX 
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APPLICATIONS  TO  WIND  TUNNEL  CONTIiACTIOJE 


In  an  ordinary  tunnel  contraction  section,  the  conditions  of  the 
linear  theory  are  not  satisfied  amruhere  in  the  contraction.  Viscous  dis- 
sipation is  important  and  there  is  a continuous  readjustment  of  the  turbu- 
lent energy  in  the  contraction  since  anisotropic  flews,  left  to  their  own 


devices,  tend  to  become  isotropic. 


Ribner  and  Tucker  hare  compared  the  results  of  their  analysis  with 
available  experiments  by  accounting  for  viscous  dissipation  in  an  approxi- 


~ A-  ~ — 


met  uo  brevy  • 


he  changes  li- 


the con traction  are  considered  to  be  a step  by 


step  process,  where  for  each  step  there  is  a change  of  energy  associajted 
with  the  stream  deformation  and  a change  of  energy  due  to  viscous  dissi- 
pation- In  most  of  the  experiments,  turbulent  intensity  measurements  ware 
made  in  the  settling  section  and  in  the  test  section.  Surprisingly  enough, 
the  experiments  are  in  reasonably  good  agreement  wi  th  the  cal  dilations - 
However,  HacFhail“'  made  measurements  on  trie  axis  all  along  the  contraction. 
While  the  ratio  of  the  lateral  turbulent,  intensity  at  the  test  section  to 
that  in  the  settling  section  is  not  too  far  from  the  prediction,  the  lateral 
velocity  ratios  along  the  contraction  are  much  larger  than  predicted,  and 
are,  in  fact,  larger-  than  would  be  obtained  neglecting  viscous  decay, 
KacPhaii  also  made  some  measurements  in  a channel  of  approximately  constant 
velocity  where  the  cross-section  was  gradually  changed  fran  a high  narrow 
rectangle  to  a low  wide  rectangle  of  the  same  area  - Again  the  development 
of  the  turbulent  velocity  components  does  not  seem  to  be  predictable  on  the 
basis  of  fluid  element  deformation.  In  addition  to  unpredicted  variations 
of  the  turbulence  intensity  in  the  deforming  section,  there  suems  to  be  an 
oscillatory  exchange  of  energy  between  the  fluctuation  components  in  the 
fixed  section  following  the  deformation . MacFhail  suggests  that  some  sort 
of  gyroscopic  action  may  be  involved  so  that  the  initial  change  of  orien- 
tation of  the  vert  icily  v=cter  is  "ot  that  predicted  by  Taylor's  linear 
theory,  in  this  concept,  the  vortex  cores  exhibit  some  of  the  character- 
istics of  30lid  bodies. 
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Sose  cossssnts  can  be  on  Haerheil^??  9Ugg*«  ti.cn  ih« i » gyroscopic 

action  may  be  present.  Tor  Solid  budioC y m.  prscsesicnal  motion  of  a 
roiatir^r  body  is  maintained  by  an  external  couple.  since  the  processional 
motion  ir.roi.vos  a time  rav-e  of  change  of  .ngular  moment u«  of  the  rotating 
body.  Tho  nean  notion  in  a wind  tunnel  contraction  (negjecting  the  wall 
boundary  layer)  is  represented  by  an  irrotational  solution.  This  ae«ina 
that  the  fluid  element  deformations  imparted  by  the  contraction  do  not 
change  the  angular  eessnius:  of  the  fluid  elements.  Therefore,  iii«  fluid 
in  the  vortex  core  should  not  process  at  least  in  the  conventional  gyro- 
scoyic  sense.  If  KacPhail 1 s measurements  are  correct  and  the  vortices 
turn  in  a way  net  predicted  by  the  fluid  deformation  theory,  then  possi- 
bly an  explanation  of  the  measurements  should  require  consideration  of 
secondary  flow  in  the  vortax  core*  as  they  p«3s  through  tho  contraction. 
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